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THE VALUE OF QUANTITATIVE INTERPRETATION 
OF GRAVITY DATA* 


D. C. SKEELS} 


ABSTRACT 


Although there is no unique interpretation of a given set of gravity data, there are 
many cases in which quantitative interpretation is decidedly worthwhile. This is es- 
pecially true in cases where the gravity data are supplemented by a certain amount of 
geological data, or where the gravity anomaly is of such a shape that the range of pos- 
sible solutions can be rather closely limited. Three examples are given of interpretations 
of actual data. 


It is a fundamental weakness of the gravitational method of geo- 
physical prospecting that the data cannot be interpreted uniquely; 
that is, for any given set of data, there is a variety of subsurface 
distributions of mass, each of which will satisfy the given data. This 
unfortunate, but inescapable, limitation has led many to the view 
that gravity data should be regarded as merely qualitative evidence 
of structure, and that quantitative interpretation is valueless and 
even misleading. 

While this view may often be justified, and while there have un- 
doubtedly been attempts to derive more from the data than they 
‘are capable of yielding, there are some cases, as I shall attempt to 
illustrate, in which the combination of gravity data with a limited 
amount of other information enables one to rule out many of the 
theoretically possible solutions, and to find one or more solutions 
which, although not unique, are consistent with all of the data and 
with geological experience, and may be regarded as good working 
hypotheses. Such an hypothetical solution need not be misleading 
if it is clearly labeled as such; and it may often prove useful in the 
extrapolation of incomplete geological data. To those who abhor ex- 
trapolation and hypothesis, it should be pointed out that many of. 


* Presented at the Denver meeting, April, 1942. 
t Standard Oil Company of New Jersey. 
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the decisions of the petroleum geologist must be based upon extra- 
polation from rather meager data. It is true that the use of the seis- 
mograph as removed some of the necessity for ‘“‘guessing,”’ but there 
are many cases in which seismic data are not available, and in which 
the geologist must base his decision on gravity data plus a certain 
amount of geological information. His decision, which may involve 
expenditures of large sums of money by his company, is usually 
based upon some hypothesis about subsurface conditions. Quanti- 
tative interpretation of the gravity data may indicate whether or 
not a given hypothesis is consistent with the gravity data; and it 
may suggest other hypotheses which fit the gravity data better, and 
the geological data equally well. The choice between hypotheses 
which are equally favored by the geophysical data should, I believe, 
rest with the geologist. 

The examples given here were taken from the files of the Standard 
Oil Company of New Jersey, and were chosen to illustrate certain 
types of information that may be derived by a quantitative study of 
the gravity data. The anomalies chosen were all of sufficient elonga- 
tion that they could be treated as two-dimensional without intro- 
ducing serious error, and this greatly simplified the calculations. For 
obvious reasons, no geographical location or other identification 
can be given for the data. 


Example I (Fig. 1) 


The available geological data are obtained from the wells shown 
on the left side of the profile. In addition, the base of the bed “C” 
was cetermined by seimsograph on the right-hand side. The contacts 
between ‘“‘A” and “B” and between “C” and “D” are unconformities. 
The outside well, at the extreme left of the profile, was dry. One of the 
wells, that on the extreme right, penetrated the granite basement. 

The gravity data show a high to the right of the field. It is ob- 
vious that the observed gravity must be due to a mass deeper than 
the seismograph horizon, and considerably steeper. The few density 
determinations which are available (from cores) indicate that the 
only important density contrast is that between the granite and the 
sediments. From geological considerations, it seems probable that 
the surface of the granite must be at least as steep as the top of the 
bed “D,” on the left flank. If we draw the granite surface conformable 
to the top of ““D,” in the vicinity of the wells, we find that the gravity 
data will be satisfied with a density contrast of 0.15, and that the gran- 
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ite on the left flank would have to have the configuration shown by 
the broken line. Combining this interpretation with the seismic data, 
we are led to predict that the producing horizon at the top of bed 
“D” should be encountered on the left flank of the structure in ap- 
proximately the position indicated. 
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Fic. 1. Use of gravity for extrapolation of geological data. 


The assumptions involved in this interpretation are: 


1. That the gravity anomaly is all due to the granite-sedimentary 
contact. A small part of the anomaly may be due to the contact 
C-D, but this cannot be very important, since, on the evidence of 
the seismograph, this surface is relatively flat on the right-hand side, 
whereas the gravity anomaly is steep on both sides. Another pos- 
sibility is that the anomaly is due to an intrusive mass within the 
basement, but, if this were true, it would be hard to explain the rather 
close correlation between gravity and structure on the left side of the 
profile. 

2. The assumption is made that the granite surface has about 
the same slope as the top of “D.’’ The granite is not likely to be less 
steep than this horizon; it may however, be steeper, in which case 
the density contrast would have to be less than 0.15, and the granite 
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structure larger in amplitude. Against this alternative, we have the 
fact that our assumed density differential of 0.15 is actually some- 
what lower than that indicated by the measured core densities; if we 
make the granite surface steeper, the discrepancy will be still greater. 

It will be noted that although “‘C” rests on “‘D,”’ on the left flank, 
our interpretation demands additional beds between ‘‘C” and “D” 
on the right flank. This is consistent with the paleontological evidence 
for there is a fairly large time break between ‘‘C” and ‘‘D.”’ 

This example is given to illustrate the possibility of using gravity 
data to extrapolate incomplete geological data. The interpretation is 
obviously not the only one possible, but is regarded as a plausible 
hypothesis, which is consistent with all of the known data. 


Example IT (Figs. 2 ana 3) 
Where geological structure is well delineated by seismograph and 
drilling, gravity data can add little to our knowledge of the structure. 
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Fic. 2. Geological, gravitational, and seismic data for a producing structure. 





A gravity survey over such a known structure may, however, give 
information about differential densities, which information will be 
useful in interpreting other gravity anomalies in the same or ad- 
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jacent areas. In this example we have the residual gravity anomaly, 
after removing the regional effect, and the geological structure as 
revealed by the seismograph and drill (Fig. 2). It will be seen that 
the gravity anomaly agrees in position with the structure, and there 
can be little doubt that the anomaly is chiefly due to density con- 
trasts within the sedimentary section, which was only partially pene- 
trated by the deepest well shown. 

One difficulty is that we do not know at which horizon the density 
contrast occurs; and, indeed, a consideration of the lithology suggests 
that there probably is no single horizon of large contrast, but that 
more probably the density increases with depth either in a continuous 
or in a step-like manner. We can, however, calculate what the density 
contrast would be if it were all located at one horizon, and if we do 
this for different horizons, the total density differential will probably 
lie somewhere between the extreme values. In Fig. 3-A the contrast 
is assumed to be at ‘“‘B,” the seismic horizon. The effect of the excess 
mass between this horizon and the sloping base line was calculated, 
and is shown as a broken line; the observed profile (residual gravity) 
is shown by the solid line. The base line was made to slope, because 
the seismograph indicates regional dip, and the effect of this regional 
was eliminated from the residual gravity profile. 

Figs. 3-B and 3-C show the calculated effects when the density 
contrast is taken to be at horizons “‘A” and “D,” respectively. For 
each horizon, that density was found which would make the maximum 
value of the calculated anomaly agree with that of the observed anom- 
aly. Horizons “D” and ‘“‘A” were extrapolated so as to be conform- 
able to “B,” which is known by seismograph. 

Finally, in Fig. 3-D the combined effect of all three horizons, with 
the density contrasts .o5, .o4, and .o4 at horizons “A,” “B,” and 
“D” is shown. 

It is obviously impossible to determine from the gravity data 
how the densities are distributed. The assumption that all of the 
density contrast is located at the deep horizon ““D” (Fig. 3-C), can 
probably be ruled out, for even if we mal< the structure considerably 
sharper on this horizon, the calculated 2:omaly will be too wide to 
match the observed profile. Fig. 3-D gives the best fit of all; but we 
could have obtained an even better fit in 3-B if we had changed the 
shape of the bed “‘A” slightly. No matter how the densities are dis- 
tributed, however, the total density contrast cannot be more than 
about 0.15, if the gravity and seismic data are both right. This is a 
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rather low figure, considering the normal increase of density with 
depth in sedimentary beds. A possibility that we have not considered 
is that there is a reversal of density contrast at some horizon, with 
an increase at another horizon. For example, if the material above 
horizon ‘‘A” were denser than that between ‘‘A” and “B” by, say 
0.10, then we might have a density contrast of about 0.25 on “B”’; 
but in this case the total density differential would still be about 0.15, 
the upper differential being regarded as negative. 


Example III (Fig. 4) 


In the previous example the depth of the disturbing mass could 
be limited to some extent, although we still had a rather wide range 
of possible depths. It sometimes happens that the shape of the anom- 
aly is such as to limit the depth rather closely. Fig. 4 shows a gravity 
anomaly which, from its shape, must be due to a steep-sided structure 
of considerable width, either an intrusion or a horst. The latter seems 
more probable from a consideration of the regional geological set- 
ting. No geological data are available in the immediate area of the 
anomaly, but the geologic section exposed in adjacent areas consists 
of a thick limestone series overlain by shales, and it is assumed that 
the only important density contrast is at the surface of the limestone. 
If this assumption be correct, we can use the gravity data to get 
some idea of the depth of the limestone on the axis of the anomaly. 
Fig. 4 shows four assumed structures, and the calculated gravity 
effect of each, compared with the observed effect. Fig. 4-A gives a 
very poor fit, 4-B fits somewhat better; Figs. 4-C and 4-D are much 
better than either of the other two. From these figures we can con- 
clude, I believe, that the limestone can hardly be deeper than 500 


‘meters on the axis of the anomaly, and is probably considerably 


shallower than this figure. 

The validity of this estimate depends, of course, upon the validity 
of the data, and of our assumption that the limestone surface is the 
only important density contrast. As far as the data are concerned, 
it should be stated that this is only one of several profiles across the 
anomaly, all of which show about the same character, so the sharp 
“corners” which reveal the shallow nature of the disturbing mass 
can hardly be accidental. As for the assumption that the limestone is 
the only important density contrast, we would feel much better if 
we had actual data on the specific gravities of the rocks exposed in 
adjacent areas. It is just possible that there is a density contrast 
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somewhere within the shale section without any corresponding change 
in visible lithology. 

The examples that have been cited are all somewhat exceptional, 
in that the anomalies are greatly elongated in one direction, and could 
be treated as two-dimensional. Where this condition is not fulfilled, 
the calculations of the gravitational effects of assumed structures are 
rather tedious, although still quite possible. Whether or not the re- 
suits will justify the time spent will depend, of course, on the problem, 
and on the amount of other data available. 

I am indebted to the officials of the Standard Oil Company of 
New Jersey for permission to publish this paper. 





















































A MECHANICAL INTEGRATOR FOR THE COM- 
PUTATION OF GRAVITY ANOMALIES* 


ARNOLD J. F. SIEGERTt 





ABSTRACT 


With the instrument here described the gravity anomaly caused by an arbitrarily 
given finite body can be computed. The integrator is operated by tracing the contour 
lines representing the body in a manner similar to the operation of a planimeter. The 
instrument is useful mainly in the indirect method of interpretation of observed gravity 
anomalies, in the computation of terrain corrections, and in the evaluation of magnetic 
anomalies. Its design is based on a formula which expresses the gravity anomaly caused 
by an infinitely thin slab, bounded by a closed curve, as a line integral over this curve. 


The gravity integrator described here was designed mainly as an 
aid in the indirect interpretation of gravity anomalies. In this method 
a geologic body is assumed, the gravity anomaly caused by it is com- 
puted and compared with the observed anomaly. If there is no agree- 
ment between computed and observed anomaly, the shape, depth, 
and density of the hypothetical body have to be changed. If there 
is reasonable agreement, the assumed body can be the cause of the 
gravity anomaly. The main disadvantage of this method lies in the 
excessive amount of computation work it requires. Even if the shape 
of the hypothetical body is oversimplified, so that its gravity anomaly 
can be expressed in terms of simple functions, the amount of work 
involved in the numerical evaluation of the theoretical gravity anom- 
aly is still prohibitive in most cases. 

The mechanical integrator here described eliminates both the 
mathematical and the numerical work, and replaces them by an 
operation as simple as the operation of a planimeter.! The instrument 
is shown in the photograph of Fig. 1, and, schematically, in Fig. 2. 
It comprises a main arm PS to which a stylus and a perpendicular 
extension arm ST are attached rigidly in the point S. The depth is 
set on the extension arm; the length of ST represents the depth. An 
auxiliary arm PT is pivoted to the extension arm in 7. Both main 


* Presented at the annual meeting in Denver, Colorado, April, 1942. 

T Stanolind Oil and Gas Company, Tulsa, Oklahoma. 

1 An integrator for the two-dimensional problem was constructed by G. A. Gam- 
burzeff, Gerl. Beitr. 24, 83, 1929. An instrument for the evaluation of certain integrals 
occurring in gravity interpretation has been described by F. Kaselitz, Zeit. Geophys. 
8 (3/4), 191, 1932. The latter can be used for the three-dimensional problem but re- 
quires plotting of intermediate results and an additional integration. (See also Geo- 
physical Exploration, C. A. Heiland, pp. 157, 269.) 
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arm and auxiliary arm slide and turn in the pivot point P. The wheel 
arm PW is perpendicular to PS and of constant length. The plane 
of the wheel W is made to form an angle a= ZSPT with PW and, 
therefore, an angle go° —a with PS. In the actual construction it turned 
out to be advantageous to attach the wheel not to the wheel arm it- 





Fic. 1. Photograph of the integrator. 


self but to a pantograph arrangement, which holds the wheel at the 
correct place and under the angle required. 

To find the gravity anomaly of a given body a contour map of 
this body is drawn on a convenient scale. The pivot point of the 
instrument is set on the point of observation, that is, the point for 








Fic. 2. Schematic diagram of the integrator. 


which the gravity anomaly is to be computed. It is held in place by a 
needle protruding from the bottom of the pivot. The stylus S is 
placed on an arbitrary point of the innermost contour line. The wheel 
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W is turned to make the wheel scale read zero, or its reading is noted 
and later subtracted from the result. The depth represented by the 
contour line is set on the depth scale ST. The depth scale and lateral 
scale must be the same; if, for instance, one inch on the map repre- 
sents 4000 ft., the depth scale has to be set to one inch for the 4000 ft. 
contour, to 3 inch for the 2000 ft. contour, etc. The stylus S is then led 
once around the contour line. The direction in which the stylus has 
to be led depends on whether the wheel scale shows increasing read- 
ings if the stylus is led around P in clockwise or in counterclockwise 
direction. It shall be assumed that the wheel reading increases if S 
is led around P in clockwise direction. To find the gravity value, 
all contour lines have to be traced in counterclockwise direction. 
It is easily seen that then the wheel reading increases when a curve 
not containing P is encircled, and decreases if the curve contains P. 
After the innermost contour line has been traced the stylus is moved 
to the starting point on the next contour line and the depth of this 
contour line is set on the depth scale. While the starting point on the 
first contour line is arbitrary all other starting points have to lie on 
a straight line through P and the first starting point and the stylus 
has to be moved on this straight line while changing contour lines, in 
order not to disturb the wheel reading. Care must also be taken not to 
move the stylus while changing the depth setting, so that the wheel 
will not turn in this procedure. After setting the depth scale, the con- 
tour line is traced, the stylus moved to the next contour line, and the 
procedure is repeated for all contour lines of the hypothetical body. 
After all contour lines have thus been traced, the wheel scale is read. 
If none of the contour lines encircled the point P, the reading J of 
the wheel is multiplied by ypt/c to obtain the value of the gravity 
anomaly at P. p is the density excess of the body over its surroundings; 
t is the contour interval; c is the calibration constant of the instru- 
ment and y is Newton’s constant. The numerical value of y, if the 
result is desired in gravity units (1/10 milligal) is: 


66.6 if¢isinkm 
107.2 if ¢ is in miles 
0.0203 if ¢ is in feet 
2.03 if ¢is ‘n 100 feet 


The above formula has to be changed slightly if some of the con- 
tour lines encircle the point P. If m is the number of these contour 
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lines, the value 2mnc has to be added to the wheel reading before 
multiplication with the above constant. The formula becomes: 


vot 
g = — (I + 2nnc). 


C 

The calibration constant of the instrument can be calculated, but 

it is simpler to find it empirically. To calibrate the instrument the 
stylus is led on a circle of radius 7 around the point P, with the depth 
setting at z, in clock wise direction. The wheel reading J is then given 


by 
= 2mcsina, inwhich sina = [1 + (r/z)?]-1/2. 


To obtain an accurate value of c, several combinations of r and z are 
chosen and the values of J read on the instrument. The sum of the 
corresponding values of 27c sin @ is divided into the sum of the in- 
strument readings to obtain c. In Table I a number of instrument 


TABLE I. EXAMPLE OF MEASUREMENT OF 
CALIBRATION CONSTANT 














r/2 27 sina 2mc sin a I (measured) 
5 5.620 667.0 (4, 8) 667.0 
I 4-443 527.4 (5, 5) 527-3 (7, 7) 527-6 
2 2.810 333-5 (6, 3) 334-1 (10, 5) 333-3 (14, 7) 333-4 
4 1.5237 180.9 (12, 3) 179.4 
c=118.7. (12, 3) means: measured with r=12, z=3. 


c is obtained by dividing the sum of all measured values J by the weighted sum of the 
values 2m sin a: 5.620+2X4.443+3 X 2.810+1.5237. The largest discrepancy between 
the computed values 27c sin a and the measured values of J is 1.5 divisions of the 
wheel scale (1 revolution= 100 divisions). For a body of density excess one and a con- 
tour interval of 1000 feet, this error would amount to only } gravity unit (0.025 milli- 


gal). 


readings are shown together with the values of 27 sin a and the values 
of 2mc sin a, computed with the value of c thus calculated. Circles 
for which the value of r/z is the same must give the same wheel read- 
ing. Comparison of the instrument readings and the values of 27c sin a 
shows the instrument to be accurate within 0.25 gravity unit for 
p=1 and ¢=1000 ft. 

In using the instrument it will, in general, be possible to place 
the point of observation so that the foot of the pivot does not cover 
any piece of a contour line. It is, for this purpose, advantageous to 
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choose the scale of the map as large as possible. This will also di- 
minish the error in tracing. If it is not possible to avoid such points 
of observation, the contour line should be traced as closely as pos- 
sible to the pivot, and the stylus then led around the pivot on an 
arc until it meets the contour line again as shown in Fig. 3. This 
amounts to neglecting the cross-hatched part of the slab represented 
by the contour line. To correct for this omission the cross-hatched 
part only is drawn on a larger scale and traced on this scale. The 


& 
FooT OF S 
PIVOT Ps) 
(ENLARGED) We 
roa; 
SV 





Fic. 3. Illustrating the procedure when a contour line 
passes under the foot of the pivot. 


depth setting has to be changed, of course, to correspond to the 
larger scale. The resulting instrument reading is added to the former 
one which is thus corrected for the omission of the cross-hatched part. 
Since the part of the contour line covered by the pivot can, in most 
cases, be considered a straight line, a correction table can be prepared 
beforehand by means of the instrument. A scale of angles can be pro- 
vided for on the foot of the pivot and the corrections are then read 
from the table which lists the instrument readings for various angles . 
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between the entrance and the exit of the contour line and for various 
depths. If P lies exactly on a contour line and the covered part of the 
contour line may be considered straight, the instruments reads 
cx[1+(ro/z)?|“? instead of cr, and the correction is, therefore, 
cw(1—[1+(ro/z)?]-"?), where ro is the radius of the foot of the pivot. 

The case of bodies whose density excess varies with depth can be 
treated easily by a slight change of procedure. Rather than reading 
the instrument only after all contour lines are traced, it is read after 
each contour line separately or after each group of contour lines repre- 
senting slabs with the same density excess over their surroundings. 
The instrument readings are then individually multiplied by yp;t/c, 
where p; is the density excess of the individual slab, or of the group of 
slabs in question. This procedure will be of value, for instance, for the 
calculation of the gravity anomaly caused by salt domes piercing sedi- 
ments of different densities. 

The application of the instrument to the computation of terrain 
corrections will decrease the time necessary for the computation from 
one hour per station to probably less than ten minutes, in cases 
where the terrain irregularities cannot be approximated by simple 
shapes like infinitely long inclined steps or ridges. If terrain correc- 
tions are to be computed, the instrument is operated on a topographic 
map which does not need to be very accurate. If the contour interval 
chosen is ¢, the contour lines are drawn at elevations ¢/2, 3¢/2, 5¢/2, 
etc., above and below the station elevation. These elevation values, 
regardless of sign, are used as z values. Contributions of contour lines 
above and below the station are to be added to the gravity value if 
elevation and Bouguer correction have been applied beforehand. 

In a similar manner, the instrument can also be used for the 
computation of the isostatic correction. 

The integrator is also useful for the computation of the vertical 
component of the magnetic anomaly caused by an arbitrarily given 
finite body, provided the magnetization can be considered homo- 
geneous, and caused by a known magnetic field, for instance, the 
magnetic field of the earth. The use of the instrument in this case is 
based on the Eétvés formulas which express the magnetic field in 
terms of the derivatives of the gravitational field: the vertical com- 
ponent H, of the anomaly, which usually is the only one measured, 
is given by 

KH, Og «H, dg 
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H, is counted positive if it increases the normal earth magnetic field, 
H,® and H, are the vertical and horizontal components of the 
normal earth field, counted positive, and x is the susceptibility of the 
body. The x axis is horizontal and directed towards the nearest mag- 
netic pole with x increasing in the direction of the pole. The z axis is 
vertical and z increases downward. p is arbitrary because it appears 
as a factor in g and cancels out. 

The contribution of H,“ can be obtained easily by graphical dif- 
ferentiation of the gravity map computed by means of the inte- 
grator. The increase of the gravity per unit distance in the direction 
of the magnetic pole is determined from the map and multiplied by 
kH, /yp. It is advantageous to begin this procedure in places on the 
map where this increase is largest, positive or negative, since this con- 
tribution is often so small that it may be neglected in all other places. 

The determination of the contribution of the vertical component 
of the earth field is based on replacing the differential quotient dg/dz 
by a difference quotient. The value of g is computed twice, with each 
value of z replaced by z+ 3Az and by z—}Az. The difference between 
these values for g is divided by Az, to obtain an approximation for 
dg/dz. The value of Az has to be chosen large enough to obtain the 
difference in the gravity values with sufficient accuracy but small 
enough to approximate the differential quotient well. A value equal to 
the contour interval ¢ should give satisfactory results. A less accurate 
value is obtained using z and z+-Az, instead of z+4Az and z—3Az, but 
even then the accuracy will probably still be sufficient for magnetic 
work. 

The question often arises whether an observed magnetic anomaly is 
caused by the same body which causes a gravity anomaly coinciding 
with it or in its immediate neighborhood. Under the assumption of 
homogeneous density and homogeneous magnetization, the magnetic 
anomaly may be determined directly from the observed gravity 
anomaly, without any assumptions concerning the shape of the body. 
The determination of the contribution of the horizontal component 
of magnetization obviously requires only a differentiation of the 
gravity anomaly. The contribution due to the vertical component of 
magnetization can be computed by means of the instrument. This pro- 
cedure is based on the theorem that the gravity anomaly at any 
point above the earth surface (assumed to be a plane) is equal to the 
gravity anomaly caused by the ‘‘ideal perturbing layer,” that is, an 
infinitely thin distribution of surface density proportional at any 
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point to the gravity value measured at this point. The gravity values, 
which would have been measured at an altitude Az above the sur- 
face of the earth, are thus computed by tracing all lines of equal 
gravity with the stylus of the instrument, keeping the depth setting 
constant at Az. The instrument reading, corrected in the usual way for 
contour lines encircling the point P, is multiplied by Ag/2my according 
to the Bouguer formula, where Ag is the contour interval of the gravity 
map. Knowing thus the gravity values at the surface and at an alti- 
tude Az, the magnetic anomaly at an altitude 5Az is determined. If Az 
was chosen small enough, this amounts to a knowledge of the magnetic 
anomaly at the surface of the earth. The magnetic anomaly thus ob- 
tained theoretically is then compared with the measured anomaly in 
order to decide whether the same body caused both anomalies. The 
procedure here described may sometimes be useful for gravity inter- 
pretation, too, since some authors believe that the vertical gradient 
of gravity is, in certain cases, more diagnostic than the gravity values 
themselves.’ 

In the same way, the gravity values at an altitude of the order of 
one mile can be computed. By subtracting the computed values from 
the measured values the influence of deep-seated structures.is greatly 
reduced. This procedure may be used instead of subtracting regional 
gradients to obtain a map on which the effect of shallow localized 
structures stands out better. This method is, of course, of qualitative 
value only, but may serve its purpose better than the subtraction of 
rather arbitrarily chosen regional gravity features, in cases where the 
regional gradient is not outstanding and simple. 


ANALYTIC BASIS OF THE DESIGN OF THE INTEGRATOR 


The fundamental operation of the instrument is the computation 
of the gravity anomaly of an infinitely thin plane horizontal layer 
bounded by an arbitrarily given curve. The value thus obtained ap- 
proximates the gravity anomaly due to a slab of finite thickness 
bounded by the same curve, within the limits of accuracy desired for 
gravity interpretation. Any given finite body can be composed of such 
slabs, with any accuracy desired. The contour lines of the body are 
the peripheries of these slabs. The accuracy obtained is very good as 
long as the thickness of the individual slabs is smaller than the 
shortest distance from the point of observation to the body. In case of 


2 Evjen, H. M. The Place of the Vertical Gradient in Gravitational Interpretations. 
Geopuysics, Vol. I, No. 1, pp. 127-136 (Jan., 1936). 
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| doubt, the computation may be repeated with a narrower contour 
: interval. For bodies with large vertical extension and steep flanks, 
another design for the instrument will be described later. 

The instrument was designed on the basis of a formula which 
expresses the gravity of an infinitely thin slab bounded by the curve 
C asa line integral over this boundary: 
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z is the depth the slab, r=r(¢) is the equation of the curve C in polar 
coordinates around the point P. The instrument performs the integra- 
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Fic. 4. Top and side view of special form of element used to derive the contour integral 
which forms the analytic basis of the design of the integrator. 
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tion if the curve C is traced with the stylus. A simple derivation of the 
principle of the instrument is based on the formula of an infinitely 
thin element of the shape shown in Fig. 4 in top view and side view. 
The element is part of a plane at depth z below the surface of the 
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earth. It is bounded by two cencentric circles and two straight lines 
intersecting in the common center of the circles. The gravity anomaly 
g due to this element at a point on the earth surface vertically above 
the center of the circles is given by 





Z Z 
= vis as a 

where y is Newton’s constant, p the density excess of the body over 
its surroundings, 7; and 72 are the radii of the circles and Ag is the angle 
between the two straight lines. If a; and a2 are the angles shown in the 
vertical section in the lower part of Fig. 4, this formula can also be 
waren g = ypt¢(sin a; — sin ae). 
This formula is the limiting case for small thickness of formula 7-410, 
p- 149, in Geophysical Exploration by C. A. Heiland. 

If the stylus S of the instrument moves on a straight line through P 
the wheel does not move. If it is led on a circle of radius r around P 
through an angle A@¢ the center of the wheel moves on a circle of 
radius PW through the same angle. The wheel partly slides, partly 
rolls. The rolling part of the motion which is read on the wheel scale 
will amount to c sinaA@. The constant c is determined by the length 
of the wheel arm PW, the radius of the wheel and the scale. If the 
stylus is led in a counterclockwise direction around the element, 
shown in Fig. 4, the wheel scale will read J=cA@ (sin ai—sin az). 
Comparing this value with the formula for the gravity g of this ele- 
ment, it is seen that 

g,= Bed i 
C 

This relation between wheel reading and gravity will now be 
shown to be true in the general case of an infinitely thin layer bounded 
by an arbitrarily given curve C if P lies outside of C. The case of P in 
C will be treated later. In Fig. 5 the area in C is shown roughly ap- 
proximated by elements of the shape considered before. To find the 
gravity at the point P, the outlines of these elements can be traced 
as shown by the arrows, the resulting instrument readings added and 
the above formula applied. It is seen, however, that the lines in the 
interior would be traced twice in opposite directions, and therefore 
may be omitted. Only the outside boundary need be traced. By in- 
creasing the number of elements the outline can be made to approxi- 
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mate the curve C more and more closely. The contributions of the 
dividing lines inside C always cancel, and to obtain the gravity due to 
the area bounded by C, therefore, ihe curve C has to be traced, and 
the equation g=Jypt/c holds for the arbitrary curve C, since it holds 
for the simple element for which it was proven initially. 








Fic. 5. Illustrating the argument used to extend the contour integral to the case of a 
thin element with an arbitrary boundary. 


If P lies inside C, the curve C is approximated by a curve C’, as 
shown in Fig. 6. C’ is identical with C, except for the piece AA’, 
which is replaced by a straight line leading into C, a circle around P 
and a straight line leading out of C, which is parallel to the first one. 
In the limit when the two straight lines are made to approach each 
other and the circle is contracted to zero radius, the area bounded by 
C’, causes the same gravity anomaly as C, but C’ does not encircle P. 
The gravity anomaly due to the area bounded by C may therefore be 
obtained by tracing C’ and applying the above formula. The two 
straight lines will not contribute to the instrument reading, because 
they are traced in opposite direction. The infinitesimal circle around P 
contributes the amount 2c. Instead of tracing C’, therefore, the 
original curve C may be traced and the value 2mc added to the instru- 
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ment reading. The value of the gravity anomaly due to an area 
bounded by C, therefore, is equal to 


yet ? — 
g = — (I — 2zc) if P lies inside C. 
C 


If m of the contour lines representing the body encircle the point P, 
is given b 
g is given by ~ 
g = — (I — 2nnc) 
C 
where 7 is the instrument reading after all of the contour lines have 


been traced. 


> 


Fic. 6. Contour used in the integration when the point of 
observation lies within one of the contour lines. 


For bodies with large vertical extent and steep flanks, the present 
instrument would necessitate the tracing of a large number of contour 
lines of nearly the same shape. A similar instrument has, therefore, 
been designed for those cases. 

The gravity value for a vertical cylinder, extending from a depth 
2, toa depth z2 and bounded by the curve C, is obtained by integration 
of the above formulas for the infinitely thin cylinder: 


7 z 
s= =j yedz Salt = wh (VP Fa VP Fad 
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To compute the integral §/r?+-22d@ the wheel, instead of being 
attached to a wheel arm PW is placed at the point 7. The wheel plane 
is kept perpendicular to PT. The instrument is otherwise constructed 
as shown in Fig. 1 and 2. If the stylus moves from the point r, ¢ to the 
point r+dr, ¢-+d¢, the wheel rolls by the amount dw: 


dw = \/r? + 2 (db — da. 


If the stylus is made to trace a closed curve C, the wheel rolls the 
amount w, where 


w = $ VF FF ap - p VF FF de 


Since +/r?+2°da is independent of ¢, the second integral vanishes 
on a closed curve, and, therefore, 


w= b VF FH ad, 


The procedure to be followed in determining the gravity anomaly is 
easily derived from the above equations which express g in terms of this 
integral. 

I should like to express my gratitude to Dr. Alexander Wolf for his 
advice in the design of the integrator and for his encouraging discus- 
sions, and to Mr. L. Evans and Mr. P. Guenther for the expert design 
and construction work. The instrument was designed and built while 
I was employed by The Texas Company, Houston, Texas. I should 
like to thank the authorities of this company for their consent to the 
publication of this paper. 














LOCAL AREAL DISTRIBUTION OF VELOCITIES 
IN THE TEXAS GULF COAST* 


B. G. SWANT 


ABSTRACT 


This paper presents velocity data obtained from long reversed profile surface 
shooting on three Texas Gulf Coast projects. Velocity determinations were distributed 
areally in sufficient detail to permit construction of ‘‘Isotime” contour maps represent- 
ing velocities to a time datum. Maps are presented showing the resulting corrected 
seismic structure. Also presented are Time-Depth curves of two profiles with curves 
of nearby well velocity surveys. Two comparisons are made of velocities at common 
panes determined by the long velocity profiles and by Time-Aé analysis of shorter 
spreads. 

Data presented show that velocities determined in the manner described are quite 
variable over local areas, but can be contoured into fairly regular patterns, and that 
if applied to the ordinary seismic data profound changes in structural interpretation 
are effected; that major discontinuities in velocity stratification may be detected; that 
in the two examples described, velocities obtained from 5400 foot reversed reflection 
profiles are higher than those obtained by Time-Af analysis of 1800 foot reversed 
profiles; and that velocity anomalies determined from surface profile shooting may 
be associated with faulting in the subsurface beds. 


For several years past Continental Oil Company has been securing 
quite a large amount of velocity information from subsurface shooting 
in the regular course of seismic operations, both from long reversed 
reflection profiles and from the shorter spread Time-A? analysis. It is 
thought that some of this material may be of interest to the profession. 

In this paper are presented chiefly data secured from reversed 
spreads of about 5400 feet, although some 4500 foot spreads are in- 
cluded. These spreads were of both the reversed basement and re- 
versed surface types. Although several different techniques were used 
in assembling and computing data from the profiles, the basic prin- 
ciple in each case reverts to the general equation: 


X2? — xX? 


V2 = —____ 
ig ~ ie 


where 
X,=distance to near end of spread 
X2=distance to far end of spread 
T,=time to nearest detector 
T.=time to farthest detector 


* Presented at the 1942 Spring meeting, Denver, Colorado. 
ft Continental Oil Company, Ponca City, Oklahoma. 
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Fic. 1. Typical Time-Depth curve as determined by a o-5400 foot reversed base- 
ment velocity profile (open circles), and Time-Depth curve from a well velocity survey 
in an adjacent area, (solid circles). 
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. Fic. 2. Velocity distribution at a time datum of 1.6 seconds, as derived from 5400 
foot reversed basement spreads, League City area, Galveston Co., Texas. Contour inter- 
val is 50 feet. 


A Time-Depth curve was drawn for each profile, and depths cor- 
responding to some chosen constant time picked, plotted on maps, 
and contoured. These maps can be considered as representative of ap- 
parent velocity variation to a time datum, and are similar to the 
“Isotime’”’ maps of Hafner.' Chiefly for experimental purposes, a 
seismic structure map was usually prepared by applying the velocity 
information to dips previously computed in the normal manner. 


1 W. Hafner, “The Seismic Velocity Distribution in the Tertiary Basins of Cali- 
fornia,” Seismological Society of America Bulletin, Vol. 30, No. 4, (Oct., 1940). 
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Fic. 3. Velocity distribution at a time datum of 1.9 seconds, as derived from 5400 
foot reversed basement spreads, League City area, Galveston Co., Texas. Contour 
interval is 50 feet. 











No attempt was made to grade the individual determinations on 
the maps. For the purpose of this paper it was thought that the atti- 
tude of contours resulting from a large number of determinations dis- 
tributed over an area should be a more direct approach to the problem. 

There are of course many ways of applying velocity information 
to seismic structure maps. The method used here attempts to make 
use of areal distribution of velocity to a Time datum. Perhaps a more 
desirable method would involve use of a geologic marker or any re- 
flecting horizon on which correlations were possible. This is a more 
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Fic. 4. Uncorrected seismic contours of a phantom horizon in the zone 3000-6000 
feet, League City area, Galveston Co., Texas. Contour interval is 50 feet. 


complex procedure in the area in question, however, so the simpler 
approach was followed. The results are presented merely for their 
possible interest, and no claims are made as to the value of the 
method used. 

LEAGUE CITY AREA, GALVESTON COUNTY 


In Fig. 1 is shown a typical Time-Depth curve as determined by 
a o to 5400 foot reversed basement velocity profile. Also shown is a 
similar curve from a velocity survey of a well in an adjacent area. It 
will be noted that the depth at 1.8 seconds as shown by the velocity 
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Fic. 5. Uncorrected seismic contours of a phantom horizon in the zone 6000-10,000 
feet, League City, Galveston Co., Texas. Contour interval is 50 feet. 


profile is about 250 feet less than that shown by the well survey. This 
should not be taken, however, as a direct indication of the error in- 
volved, since the surface profile location is in a down-dip position with 
respect to the well survey. 

In Fig. 2 is shown the velocity distribution at a time datum of 1.6 
seconds. Fig. 3 is the same type of map at a datum of 1.9 seconds. The 
similarity of the two maps is quite apparent, the deeper beds showing 
the greater velocity variation. 

In Figs. 4 and 5 the seismic structure as represented by contours 
of averaged dips for zones of, respectively, 3000 to 6000 feet, and 6000 
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Fic. 6. Seismic contours of a phantom horizon in the zone 3000-6000 feet, corrected 
for the velocity distribution shown in Fig. 2. League City area, Galveston Co., Texas. 
Contour interval is 50 feet. 


to 10,000 feet, is shown. The contour interval is 50 feet. Dips were 
computed and averaged, using a constant velocity for each zone to 
convert times to feet. All shooting was of the continuous reversed 
profile type, with spread lengths usually 1800 feet. 

Fig. 6 shows the seismic contours of Fig. 4 corrected for the velocity 
distribution shown in Fig. 2. Fig. 7 is a similar map of the deeper struc- 
ture—the map of Fig. 5 corrected for the velocity distribution shown 
in Fig. 3. As to be expected, the closed velocity lows, one just south of 
Clear Lake, north central part of the map, and the other northeast of 
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F ; Fic. 7. Seismic contours of a phantom horizon in the zone 6000-10,000 feet, cor- 
rected for the velocity distribution shown in Fig. 3. League City Area, Galveston Co., 
Texas. Contour interval is 50 feet. 


the League City field, southwest part of the map, convert the seismic 
_ contours into closed structural highs on the corrected maps. The addi- 
tion of the faults in the southwest portion of the corrected maps is an 
effect of the velocity distribution on the contoured interpretation. 
This may seem to be manufacturing structures out of thin air; 
however, the subsurface map of Fig. 8, which is one of the many pos- 
sible contour interpretations of Middle Frio structure as indicated by 
the various wells, shows a faulted structure to lie under the large 
“‘manufactured”’ seismic closure northeast of the League City field. 
It may be open to question whether one should rely too heavily on 
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Fic. 8. Subsurface map of League City area, Galveston Co., Texas, contoured on the 
top of the Lobit (Middle Frio) sand. Contour interval is 50 feet. 


this apparent relationship, but nevertheless it demonstrates what 
certainly may and perhaps does happen to some degree. 

Although it can well be argued that faults not detected by the 
seismograph survey are actually responsible for the failure of the un- 
corrected dip map to check the subsequent well data, it should be 
noted that the velocity distribution as mapped is of such a nature as 
to improve the tie. Referring to Fig. 8, it is seen that on the Middle 
Frio “Lobit” sand the Conoco #1 Stewart Title well is 281 feet higher 
than the Phillips #1 Brittnacher located immediately northeast of the 
League City field. In Fig. 5, the uncorrected dip map, the relative 
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Fic. 9. Velocity distribution at a time datum of 1.9 seconds, as derived by Time- 
At analysis of 1800 foot reversed spreads, four spreads being used for each determina- 
tion. League City area, Galveston Co., Texas. Contour interval is 50 feet. Compare 
with Fig. 3. 
relief is seen to be only 100 feet between these two wells. On the cor- 
rected map of Fig. 7 the relief has been increased to 200 feet between 
the wells. Using the same attack on the relative relief between the 
Continental #1 Lobit and #1 Davis wells, it is seen that the well logs 
show 280 feet, while the corrected and uncorrected maps show 160 and 
200 feet respectively. 

Additional evidence of a deep-seated basement high here is indi- 
cated by the rapid thinning of sediments toward the bottom of the 
hole in the Continental #1 Lobit well on the upthrown side of the fault. 
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Fic. 10. Velocity distribution at a time datum of 1.3 seconds, as derived from 
5400 foot reversed basement spreads, West Manvel Area, Brazoria Co., Texas. Contour 
interval is 50 feet. 


There was little evidence from the seismic work, done at about the 
time the League City field was discovered, to indicate the fault which 
was cut by the last mentioned well, although evidence was noted of 
other possible faults shown in Figs. 5, 6 and 7. 

In comparing the subsurface map of Fig. 8 with the velocity dis- 
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Fic. 11. Velocity distribution at a time datum of 1.6 seconds, as derived from 5400 
reversed basement spreads, West Manvel area, Brazoria Co., Texas. Contour interval 
is 50 feet. 


tribution shown in Figs. 2 and 3, as previously indicated, the large 
velocity low northeast of the League City field and the prominent 
faulted structure seem to be related. Also it should be noted that the 
contours place the League City field proper on a velocity high trend, 
while Dickinson is outside the velocity control points and nothing can 
be said as to the relation there. 

Fig. 9 is a map showing the velocity distribution at a time datum 
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Fic. 12. Velocity distribution at a time datum of 1.9 seconds, as derived from 
5400 foot reversed basement spreads, West Manvel area, Brazoria Co., Texas. Contour 
interval is 100 feet. 


of 1.9 seconds as indicated by Time-Af analysis of 1800 foot reversed 
spreads, four such spreads being used for each determination. It will 
be noted that in comparing this map with that of Fig. 3, the general 
features are similar, although the prominent velocity low northeast of 
the League City field has been shifted eastward by the A¢ analysis 
determinations, and the amount of closure lessened. The slightly dif- 
ferent distribution of the determinations between the two maps may 
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Fic. 13. Velocity distribution at a time datum of 1.9 seconds, as derived by Time- 
At analysis of 1800 foot reversed profiles. West Manvel area, Brazoria Co., Texas. 
Contour interval is 100 feet. Compare with Fig. 12. 


account for some but certainly not all of this. In comparing 27 com- 
mon points on these two maps it is found that the Time-A? analysis 
velocities average 77 feet per second slower than those from the 
longer profiles. Considering the errors inherent in each determination 
and the comparatively small number of observations, this difference 
of only one per cent in overall velocity seems more or less accidental. 
However, this bears out the writer’s conclusions, after making a num- 
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Fic. 14. Uncorrected seismic contours of a phantom horizon in the zone 2500-5000 
feet, West Manvel area, Brazoria Co., Texas. Contour interval is 50 feet. 


ber of such comparisons in various areas, that, ‘‘In velocity deter- 
minations from surface shooting, the longer spreads tend to give 
higher velocities than do the shorter spreads.” : 


WEST MANVEL AREA, BRAZORIA COUNTY 


In Figs. 10, 11, and 12, the contoured velocity distribution at a 
time datum of 1.3, 1.6, and 1.9 seconds, respectively, is shown. The 
contour interval is 50 feet on the first two and 100 feet on the third. 
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Fic. 15. Uncorrected seismic contours of a phantom horizon in the zone s5000- 
10,000 feet, West Manvel area, Brazoria Co., Texas. Contour interval is 50 feet. 


The determinations were made from o to 5400 foot reversed base- 
ment spreads, as was the case in the League City area. The chief 
feature here is the prominent velocity low through the central part of 
the area, with high areas building up to the south and north. Except 
for minor details and the increased gradient with depth, the three 
maps present very much the same picture. 

Fig. 13 shows the velocity distribution at the 1.9 second datum 
of the south one-half of the area as determined by Time-A¢ analysis in 
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Fic. 16. Seismic contours of a phantom horizon in the zone 2500-5000 feet, cor- 
rected for the velocity distribution shown in Fig. 10. West Manvel area, Brazoria Co., 
Texas. Contour interval is 50 feet. 


the same manner as in Fig. 9. This was done to afford a check on the 
velocity high in the southwest part of the area. As may be seen in a 
comparison with Fig. 12, the check is very good. The sharp low in 
Section 75 may be discounted somewhat because of lack of supporting 
values. 


Again a comparison of the total average velocity represented by 
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Fic. 17. Seismic contours of a phantom horizon in the zone 5000-10,000 feet, cor- 
rected for the velocity distribution shown in Fig. 11. West Manvel area, Brazoria Co., 
Texas. Contour interval is 50 feet. 


these two types of determinations, 22 common points, shows that the 
Time-At analysis of the shorter spreads gives the slower velocity by 
182 feet per second, or about 23 per cent. 

In Figs. 14 and 15 are shown the seismic contours, uncorrected for 
velocity distribution, of a zone from 2500 to 5000 feet and a zone from 
5000 to 10,000 feet, respectively. Dips on which these contours are 











VELOCITIES IN THE TEXAS GULF COAST 385 











































































































we tN 
g dl . ie 
¥. y a} - “he ait 
0 L oe; “ SR--s, i 
“Ta8 “Tete s > ' . 
7 76 a of See. am, 
a 4 a a Fou D — z. 
Ps Cy ~ \. 
° ° ail e' ? acne 8 
Ma ae 8 olan ‘ 
en Tas ‘ . 
° ” ue 8 Lat pd i Py \ 
j ans - H 
° \ Pa ; ‘ 
° ° ’ f 
=o, 9 \ bs ? / i 
a6 io Kran yi Bee ; y 
a tT! Me _* al E a ag 
a ; 
b por, z — 
te. RS a a 
7 an 4600 
’ _ rs 
 - Tae 
<<" 4 \ \ 66 +5000 
— 
J = anca: WEST MANVEL 
EE ~ COUNTY: BRAZORIA STATE: TEXAS 
¥ 5 a, .4. 100° 
OATUM: HETEROSTEGINA 
CONTOURING INFLUENCED SOMEWHAT BY 
” Pb $ paizee SEISMOS DATA. 


























Fic. 18. Subsurface map of the West Manvel area, Brazoria Co., Texas, contoured 


on the Heterostegina, but with the contouring in the absence of well control influenced 
by the seismic data. 


based were computed in the usual manner and averaged into the two 
zones.” Both record quality and contouring convenience serve to locate 


* The sub-sea values assigned to the contours on these maps, as well as to those in 
Figs. 4, 5, 6, 7, 16, 17, and 20, are “phantom” in nature, and are not to be used as 
absolute values in comparing corrected and uncorrected maps. Neither are the values 
continuous on opposite sides of faults not completely delineated by the seismic survey. 
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the various faults shown. The Mykawa field is a short distance off the 
northeast corner of the map, and the approximate outline of produc- 
tion at Manvel is shown by the irregular dotted line in the east-central 
portion of the map. 

Figs. 16 and 17 show the seismic contours corrected for the velocity 
distribution. Although the general appearance remains the same, it is 
seen that the pronounced velocity high in the south central part of the 
area has shifted the structural high westward in Fig. 17, and also made 
unnecessary the fault through Sections 64 and 73 of Fig. 15. As an 
examination of Heterostegina points on the various dry holes in the 
area will show, neither map correctly indicates the structural relief in 
all cases. The uncorrected map is better in this respect than the other. 
No map was prepared with the seismic data restricted to the Hetero- 
stegina horizon alone. 

In Fig. 18 is one interpretation of the Heterostegina subsurface of 
the area. The contouring between well control points has been in- 
fluenced somewhat by the seismic data, especially as to the position of 
certain of the faults. Keeping in mind this fact, there is some doubt as 
to the actual relation of the apparent similarity of the velocity low 
area and the area between the two parallel northeast-southwest faults. 


The assumption does seem justified, however, that there is some real 
connection between the two, whether it be through the basic seismic 
data’ or from a structural angle. 

The contours as shown in Figs. 10, 11 and 12 suggest that both 
Mykawa and Manvel may lie in low velocity areas. Lack of data 
makes this only a suggestion, however. 


SOUTH MARKHAM AND BLESSING AREAS, MATAGORDA COUNTY 


Since this area was shot as two or three separate projects with 
accompanying changes in technique, determinations are from both o 
to 4500 foot reversed surface, and o to 5400 foot reversed basement 
spreads. The two types are not distinguished on the maps, as a careful 
examination seemed to indicate that this factor had little bearing on 
the final results. 

In Fig. 19 is shown the velocity distribution at a time datum of 
2.2 seconds. Contours are at an interval of 100 feet and, if taken 
literally, indicate a very irregular velocity pattern. Especially prom- 
inent are: (1) the long velocity low in the northern part of the area; 


3 By this is meant that the Af information from the records is used in both velocity 
and dip computations at any particular location, and therefore might result in anomal- 
ous values for both the dips and the velocities. 
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(2) the velocity high and adjacent lows near the Blessing field, west 
central part of the map; and (3) the velocity increase to the north 
away from the West Citrus Grove field. It is unfortunate that de- 
terminations were not secured to support the sharp velocity low at 
what is now the Blessing field. Unsupported, this lone point must be 
questioned, along with any tempting suggestions as to the relation of 
the low to the field. 

Fig. 20 is a map showing a contoured interpretation of the normal 
seismic dip data averaged for the 6000 to 11,000 foot zone. As in the 
Manvel area, record quality and contouring convenience account for 
the numerous faults shown. It will be noted that the gentle faulted 
structure in the John Partin Survey falls within a velocity low shown 
in Fig. 19. The same is true of the seismic high in the J. Filly and P. 
Burnett Surveys, while the structure in the J. Updike Survey is at a 
velocity high. There seems to be a possibility that other of the faults 
may be related to some of the isolated velocity lows. 

No map was prepared showing normal seismic data corrected for 
the velocity distribution in this area; however, it is obvious that the 
effects would be quiie pronounced. Closure on the gentle faulted struc- 
tures in the north part of the area, located in velocity lows, would be 
increased, and the Blessing structure shifted southward. Taken 
literally, all the closed velocity lows would tend to produce structural 
highs and the velocity highs to produce structural lows. 

Fig. 21 is a contoured interpretation of Frio subsurface points at 
the various scattered wells. No attempt was made to show the details 
of the fields. Seismic dips have been used to a great extent in the areas 
barren of well control, which makes a critical comparison with the 
seismic map of little value. It is to be noted, however, that reasonable 
contouring of the two types of data on one map is quite possible. The 
Blessing field structure seems slightly displaced to the north on the 
seismic map, possibly due to insufficient shooting, although a critical 
examination will show the displacement to be much less than it looks 
on this map- 

Fig. 22 is a Time-Depth curve of a velocity survey in a nearby 
well, with one obtained from two o to 4500 foot reversed surface pro- 
files. The outstanding feature is the abrupt velocity decrease at about 
1.6 seconds’ time on each curve. Again these two curves cannot be 
compared directly as to the velocities and depths they represent be- 
cause of the difference in location; however, the similarity of shape 
does seem significant. It is of interest to report that practically all of 
the velocity determinations made in this area give Time-Depth curves 
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Fic. 19. Velocity distribution at a time datum of 2.2 seconds, derived from 5400 
foot reversed basement spreads and 4500 foot reversed surface spreads, South Mark- 
ham and Blessing area, Matagorda Co., Texas. Contour interval is 100 feet. 
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Fic. 20. Uncorrected seismic contours of a phantom horizon in the zone 6000~ 
10,000 feet, South Markham and Blessing area, Matagorda Co., Texas. Contour inter- 
val is 50 feet. 
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Fic. 21. ‘Subsurface map of the South Markham and Blessing area, Matagorda Co., 


Texas, contoured on the Frio. Contouring greatly influenced by seismic data. Contour 
interval is 100 feet. 
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Fic. 22. Time-Depth curve as determined from two o—4500 foot reversed surface 
profiles (open circles), and Time-Depth curve obtained in a velocity survey in a nearby 
well (solid circles). South Markham and Blessing area, Matagorda Co., Texas. 
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of this same characteristic shape. The additional inflection point at 
about 2.6 seconds is also characteristic of other curves in the area, 
although less significance is attached to it because of the wider scatter- 
ing of points at the greater depth and lack of a definite subsurface tie. 

The log of the well in which this velocity survey was made shows 
the inflection point at about 6000 feet to be the top of the Marine 
Oligocene. The fact that this section is here mostly shales, while the 
Miocene immediately above is chiefly sands, is thought to be the 
proper explanation of the inflection in the curve. 

This phenomenon was noticed and properly interpreted in the 
curves from the surface velocity profile shooting by the field crew in 
routine procedures, and without previous knowledge of the well survey 
data. This seems significant, especially since the usual tendency is to 
smooth out the curves as much as possible and to discount points 
that do not fit as being due to errors inherent in the method. 


SUMMARY 


1. Data from three Texas Gulf Coast projects show that velocities 
determined from long reversed reflection profiles and from Time-Aé 
analyses of shorter profiles are quite variable over local areas, but can 
be contoured into fairly regular patterns. If applied to the ordinary 
seismic data profound changes in structural interpretation are effected. 
The evidence presented does not indicate the validity of such a pro- 
cedure, nor is it necessarily condemned. 

2. Velocities determined from long reversed reflection profiles in 
Matagorda County, when plotted as Time-Depth curves or as some 
other function of time or depth, may indicate clearly major discon- 
tinuities in velocity stratification. In the example presented, one such 
discontinuity is definitely identified as the Mzocene-Oligocene contact. 

3. Intwoareas it is shown that velocities obtained from 5400 foot 
reversed reflection profiles are higher than those obtained by Time-Aé 
analysis of 1800 foot reversed profiles. There is noted, however, a 
general agreement in the contoured areal distribution. 

4. Velocity anomalies determined from surface profile shooting 
may be associated with faulting in the subsurface beds. 

The author wishes to acknowledge his gratitude to the Conti- 
nental Oil Company for permission to publish this paper, and his 
appreciation of the original work of C. A. Putnam and R. W. Carson 
in securing the field data. 














AN ANALYTICAL METHOD OF MAKING 
WEATHERING CORRECTIONS* 


M. B. DOBRINt 


ABSTRACT 


A method of weathering is described by which intercept times can be rapidly and 
accurately computed from first arrival times without the plotting of time-distance 
curves. The velocities are determined by a mechanical procedure, based on least squares 
theory, which normally requires no exercise of judgment on the part of the computer. 
The application of the method to actual field set-ups is illustrated by sample calcula- 
tions. 


INTRODUCTION 


In the conventional method of determining seismic weathering cor- 
rections from intercepts on first arrival time curves, it is necessary to 
ascertain the velocity of seismic waves in the first high-speed zone 
beneath the weathered layer. This is generally done by plotting the 
first arrival times against the shot-detector distance and estimating 
the speed from the slopes of the curves thus obtained. The usual pro- 
cedure for determining these slopes is graphical and consists in averag- 
ing the deviations by eye of each plotted point from the best line that 
could be drawn through them all. 

This subjective process is often difficult and time-consuming, since 
in most cases it involves the use of trial-and-error methods. A com- 
puter, for this reason, must have had considerable experience before he 
can make reliable weathering corrections for records shot in country 
with any but the most straightforward surface conditions. 

The use of standard least square methods for determining the best 
line through a series of points would remove this undesirable subjec- 
tive factor from the weathering calculation, but these are generally so 
laborious that they would not be practical for routine work. However, 
by utilizing certain characteristics of ordinary reflection or refraction 
set-ups, such as the regular spacing of detectors and the coincidence of 
detector positions for shots from opposite directions, a least squares 
method can be set up which enables the whole weathering calculation 
to be performed mechanically and which is also more rapid and accu- 
rate than the ordinary procedure. By this method it is possible to 
calculate the intercepts directly from the first arrival times and the 


* Presented at the Spring 1942 meeting, Denver, Colorado. 
} Gulf Research and Development Company, Pittsburgh. Now at the Naval 
Ordnance Laboratory, Washington, D. C. 
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distances without the necessity of plotting the arrivals or of determin- 
ing the velocity of the high-speed zone explicitly. 


THEORY 


The present method is based upon the principle that the relative 
values of the intercept times as determined by the shooting in opposite 
directions across a spread should be independent of the direction of 
shooting. If there is no weathering or equal weathering under the shot 
points at opposite ends of a reverse spread the intercept times at each 
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TIME-DISTANCE PLOT OF TYPICAL SPREAD 


Fic. 1. Time-distance diagram of a typical spread. 


detector position should be identical for the two shots. If there is 
unequal weathering, the two intercept times at each detector position 
should differ by constant amounts. 

In the analytical procedure which forms the basis for this method, 
the velocities are determined by performing simple arithmetical opera- 
tions, such as multiplying by small digits, upon the differences, at 
each detector, between the arrival times associated with the two 
shots. These operations effectively weight the various differential 
times so that each makes its proper contribution to the velocity. This 
quantity is not expressed directly but instead is put into the form of a 
time parameter, s/V, where s is the separation of the shot points at 
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opposite ends of the spread. From this parameter and the sum of the 
arrival times at each detector position, the average intercept time is 
computed by simple subtraction. 

Fig. 1 shows a time-distance plot for a set-up of m detectors, 
spaced a constant distance / apart, located between the shot points A 
and B, which are separated by the distance s. Detector 1 is a distance 
a from A. The arrival times as recorded at the detectors 1, 2,...” 
are to be designated as T41, Ta2,... Tan for shots from A and as 
Tz, T2,...- Ten for shots from B. The actual velocity, V, can be 
obtained by averaging the apparent velocities determined by the 
arrival times from each direction; and this can be best represented by 
two straight lines of slopes 1/V (representing the shot from A) and 
—1/V (representing the one from B). These are the ideal time 
distance curves. Where the ideal curve for a shot from A crosses the 
point midway between detector 1 and detector , we will designate the 
time as J». Where the line from B crosses the midpoint we will desig- 
nate the time as Ty’. It is of course not necessary that this midpoint 
be equidistant from A and B. 

Now, if at the kth detector the arrival times were to lie on these 
ideal lines we would have 


ed (ee ey ee 
Ak = 5 [4 — ( - | 
= T) — (n+ 1 — 2k)l/2V 


theoretically. But the actual value of 74; will differ from this by an 
amount 


(x) 


€Ak = Tak — To -++ (n + = 2k)l/2V. (2) 
Similarly, shooting from the opposite direction 
To = Ty abo — (k- ot] 
V 2 (3) 


= T) + (n+ 1 — 2k)l/2V 
theoretically, and 
esr = Tar — To — (n +1 — 2hk)l/2V. (4) 


To determine the best velocity, V, from the plotted points we 
make use of least squares theory and choose V as that velocity which 
will make the sum of the squares of the e values a minimum. That is, 
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>> & = aminimum, or (5) 
d(e?) de 
fe age oe ee (6) 


Applying (6) to (2) and (4), 


o= 2 { Tar—Tot+(n+1—2k)l/2V } (n+1—2k)l/(—V?2) 
(7) 
+ > { Tex—To'+(n+1—2k)l/2V } (— [n+1—2k]l/(—V?). 


The To and Jo’ terms drop out because for each positive value 
of (n+-1—2k) there is a negative value that cancels it, so that 


o= - (Tar)(n + 1 — 2k) + [(n +1 — 2k)1]?/2V 
mn (8) 
+ (= Tar)(n +1 — 2k + [(n +1 — 2k)]?/2V. 


Solving this for V, 
v= — dL (n+ Sow ane [|S (Tax — Tai)(n +1 — 2k)l. (9) 
1 1 


This furnishes the velocity at the base of the weathered layer. 
To determine the intercept time at the Ath detector for the arrival 
from shot point A one sets 


bax = Tat — {a + (~- 1)1]/V. (10) 
Similarly, 
be. = Tar — {s— [a + (k — 1)l]}/V. (11) 


Now the best intercept time for the kth detector location will be the 
average of b4, and Ddzx: 


b= (barter) /2= (Tar tT px)/2—s/2V (12) 


| - (T44—Tpx)(n+1— 2k) 
§ 1 
=4(Tai+ Ta) — | —— - (13) 


. > (n+1—2k)? | 
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APPLICATION 


Eq. 13 enables one to compute the intercept time for a detector 
location in terms of the arrival times, the number of detectors, the 
distance between adjacent shots, and the detector separation. In 
routine shooting, all of these quantities are generally fixed except for 
k, the ordinal detector number, and the arrival times, T4, and Tax. 

Actually, the (n-+-1—2k) terms are weighting factors which de- 
termine the relative contributions of the arrivals at each detector in 
the determination of the velocity. When k approaches zero or n, the 
absolute value of (n-++-1—2k) will be greatest, indicating that those 
arrivals nearest the shot point are weighted most heavily in this 
computation. 

Let us apply Eq. 13 to the determination of intercept times for 
the case of four detectors between the shot points. This is a common 
setup for a six detector spread, since the first detector is often used 
to record the up-hole time, and the last detector is placed at the next 
shot point. One evaluates the s/V term by a summation from k=1 
to k=n=4. 


I 
ae T1 = T 21 T a2 ae T Be 
== { 3( ) + ( ) (29) 


— (Fag > Tay) — Tie Tx) | /20l. 
Or, rearranging terms, 


I 
a = t3l(Taa— Tas) + (Tor — Tard] 


+ [(Tp2 — Tas) + (Tas — Tx) |} /20l (15) 


= [3(ATs + AT}) + (AT2 + AT3) ]/201 
so that . 
by = 2(Tax + Tax) — = {3[AT, + AT] + [AT2 + ATs]}. (16) 


A similar calculation for the case of six detectors gives 


a = [5(ATs + AT:) + 3(ATs + AT2) + (ATs +.AT3)]/70l (17) 


the intercept times being calculated as in (16) for the four-detector 
case. 
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To illustrate the application of the method, Table I is reproduced. 
It is the calculation that was used in the actual weathering determina- 
tion for a six-trace spread, with four detectors between the shot 
points. Fig. 2 shows how the weathering would be carried out with 
the same data if the conventional method, employing a time-distance 
curve, were used. The intercept values indicated for each detector 
position on Fig. 2 agree closely with those obtained in Table I. 
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Fic. 2. Time-distance plot and conventional solution. 
of weathering for the data shown in Table I. 


TABLE I. EXAMPLE OF LEAST SQUARES COMPUTATION 
OF INTERCEPTS 








Shot B/Order of 


k Shot A . . 3 Operation 





Station 49+50 52+00 54+50 57+00 59+50 62-+00 
by 50 60 63 55 47 48 = (11) 
TartTp —s/V 121 127 III 94 (10) 
Tart Tor 835 841 825 808 
Tak 59 205 35t 485 620 762 
T pr 764 630 490 340 188 50 
AT; 425 139 145 432 
AT,+AT, 857 
3(AT1+AT%) 2571 
AT2+AT3 284 
3(AT1 +474) +(AT2+AT3) 2855 
s/V 714 
Times are in milliseconds 





The intercept times at the two end detector positions, which are 
coincident with the two adjacent shot points, were obtained by sub- 
tracting the s/V value from each observed arrival time. It is of inter- 
est to observe that the velocity is nowhere used explicitly or even 
calculated as such. However, if the velocity is desired one merely 
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divides the s/V value into the length of the spread, s. If the distance 
between the shot on either side and the detector nearest it is changed, 
the calculation is not affected unless s changes also. 

For the more usual cases where a larger number of detectors is 
employed, the numerical coefficients will be somewhat larger and more 
numerous. Table II gives the (n-+-1—2k) coefficients for each set of 
detector positions in spreads of from four to eleven detectors. To 
illustrate its application for the case of ten detectors it should be 
noted that the outermost pair of AT values is multiplied by 9, the 
next pair toward the center by 7, etc. The last term in the table evalu- 
ates the coefficient in the numerator of equation (9). 


TABLE II. COEFFICIENTS FOR AT VALUES IN LEAST 
SQUARES WEATHERING METHOD 








k 
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As the number of detectors increases, a point is reached where 
the method becomes too cumbersome to be applied conveniently. 
The maximum number, in practice, is to be governed by the ac- 
curacy desired in the results and in the complexity of the surface 
conditions in the area being worked. 
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AN APPLICATION OF SEISMIC SURVEYING TO THE 
LOCATION OF BAUXITE IN ARKANSAS* 


J. A. GILLIN,{ LORENZ SHOCK,j{ ann E. D. ALCOCK7 


ABSTRACT 


The exploration for bauxite has been greatly stimulated by the increased demand 
for aluminum. The geology of the formation of bauxite in Arkansas is closely related 
to the syenite exposures on the old land surface at the end of Midway time. A refrac- 
tion seismic survey mapped the attitude of the high velocity formations identified as 
syenite or Paleozoic rocks. By combining this map with a map of the Midway from scat- 
tered core holes, a map of the old land surface at the end of Midway time was obtained. 
By interpreting this map, future explorations by core drills can be guided to the most 
favorable locations. 


The increased demand for aluminum in the war program, and the 
lack of shipping space in which to import bauxite from South America, 
has greatly stimulated the exploration for new bauxite deposits in 
the United States. Approximately 90% of the present bauxite pro- 
duction in the United States comes from a small region in the central 
part of the state of Arkansas.’ The deposits of bauxite in this region 
are usually of large size and compactly located in a small area. There 


are other deposits known to exist in the Southeastern states of the 
United States, which are being worked on a small scale. These de- 
posits, however, are generally of small extent and rather widely scat- 
tered, hence from the standpoint of an immediate increase in our 
reserves, the Arkansas bauxite region offers the most promise for 
exploration. 

The commercial bauxite deposits in Arkansas occur in Saline and 
Pulaski Counties near Little Rock in the central part of the state, 
(Fig. 1). The known deposits are confined to two small areas approxi- 
mately fourteen miles apart, the Sweet Home area lying about four 
miles south and southeast of Little Rock and the Saline County area 
centered around the town of Bauxite. In both these areas the deposits 
lie close to an outcrop of nepheline syenite. The syenite occurs as an 
intrusion through the Paleozoic basement rocks. The composition of 
the syenite in both outcrops is identical, which might indicate that 
they were cupolas of the same underlying syenite mass. 


* Presented at the Denver meeting, April, 1942. 

{ National Geophysical Company, Dallas, Texas: 

1 Bramlette. The Geology of the Arkansas Bauxite Region. Arkansas Geological 
Survey Information Circular No. 8. 
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Fic. 1. Arkansas bauxite region, Grant, Saline, and Pulaski Counties, Arkansas. Known 
deposits are confined to the vicinities of the towns of Sweet Home and Bauxite. 





J. A. GILLIN, LORENZ SHOCK, AND E. D. ALCOCK 


In Information Circular No. 8 of the Arkansas Geological Survey, 
Bramlette concludes from his studies that the bauxite in Arkansas 
was formed as an alteration product of nepheline syenite. This theory 
was first developed and established as a practical working principle 
by W. J. Mead? in 1912. Bramlette’s evidence also indicates that the 
surficial alteration took place during a single geological period, the 
transition period between the Midway and the Wilcox, hence, the 
bauxite deposits are confined to the old land surface existing during 
this period. Several deposits of bauxite have been found at the con- 
tact zone between the Midway and the Wilcox; however, the only 
commercial bauxite is that which is found lying on or near the nephe- 
line syenite. Granting these conclusions, the major condition for the 
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Fic. 2. Generalized cross-section of the Arkansas bauxite region, illustrating the 
occurrence of bauxite on or adjacent to portions of the syenite mass not overlain by 
Midway sediments. 


presence of commercial bauxite at any location is that the nepheline 
syenite should have formed part of the land surface during the transi- 
tional period between the Midway and the Wilcox ages. Furthermore, 
since the bauxite was developed as an alteration product from the 
syenite through leaching by ground waters during the Midway-Wil- 
cox interval, it follows that bauxite will be found only on or immedi- 
ately adjoining syenite not covered by Midway sediments. Fig. 2 
is a generalized cross-section of the Arkansas bauxite region which 
illustrates some of these points. 

The physical characteristics of bauxite ore vary considerably and 
are quite similar to those of adjacent sediments, hence geophysical 
methods do not hold much promise for the detection of the presence 
of the bauxite itself. Furthermore, the quality of the ore can only be 
determined by chemical analysis. The only successful solution for the 
location of bauxite to date has been the use of a drill followed by a 
chemical analysis of the samples. 

A map of the old land surface during the Midway- Wilcox transi- 


2 Econ. Geol., Vol. X, pp. 28-54 (1915). 
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tional period would be of great value, however, to the exploration 
program. Such a map would outline the only areas in which bauxite 
could have been formed and would eliminate the necessity for drilling 
outside these areas. This map would also indicate areas of higher 
probability for the presence of bauxite within the areas in which 
bauxite could have been formed. 

The problem of mapping subsurface structures, or old land sur- 
faces, is particularly adapted to solution by use of the seismograph. 
Due to the shallow depth of the surface to be mapped (not exceeding 
1000 feet) the problem is more readily solved by use of the refraction 
seismograph than the reflection seismograph. 

The velocity contrast between the syenite and the overlying sedi- 
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Fic. 3. Refraction spread used in the mapping of the Syenite-Paleozoic surface. 


ments is sufficiently large to allow the mapping of the attitude of the 
syenite with the refraction seismograph. Since this surface is highly 
irregular the meagre core drill data available are useless for mapping 
the syenite, except in the immediate vicinity of the core holes. Noth- 
ing was known about the velocity contrast between the Midway and 
the Wilcox groups; however, this is a depositional surface, therefore 
nearly regular, and sufficient core hole drilling data are available to 
furnish a fairly reliable map of this surface. By combining the seismo- 
graph map of the syenite and the core drill map of the top of the Mid- 
way, a map of the old land surface existing at the end of Midway time 
can be made. 

A topographic ridge cut by several streams running between 
Sweet Home and Bauxite offers possibilities for the presence of a 
buried ridge of syenite connecting the two outcrops whose top might 
lie above the Midway zone. 

In June, 1941, a limited seismic survey of the Arkansas bauxite 
region was made by the National Geophysical Company for the 
benefit of the Reynolds Mining Corporation. The survey was con- 
ducted with a National Geophysical Company ten trace recording 
seismograph. Each trace was actuated by the output of one seismo- 
meter. The instrumental arrangement used is shown in Fig. 3. Using 
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this arrangement data were recorded in opposite directions for each 
spread. Since the computations were based on the first arrivals of 
energy from the Syenite-Paleozoic surface the distance from the shot 
point to the spread varied with the depth of this surface. The Syenite- 
Paleozoic surface is rather irregular and it is often difficult to even 
approximate the depth before recording. In some cases it was neces- 
sary to shoot at several distances, since the first distance used was 
not sufficient to give the refraction from the Syenite-Paleozoic surface 
as a first arrival. The longest distance used was 4000 feet. If the 
Syenite-Paleozoic refraction did not record as a first arrival at this 
distance it lay at a depth greater than 1000 feet and therefore was 
well below the top of the Midway. When the distances needed to be 
increased on any spread, new holes were drilled since the depth re- 
quired was only 10~20 feet and it was faster to drill new holes than 
to move the instruments. 

The computations were made using the usual refraction formulae 
for first arrival of the refracted energy. Weathering corrections were 
made using “‘up-hole times” at the shot points and short refraction 
profiles at the instrument spreads. 

The vertical velocity for the low speed layers lying above the 
Syenite-Paleozoic surface was measured by ‘‘up-hole times” in several 
deep holes drilled throughout the area which gave a fairly constant 
value of 5600 feet/sec. In addition, the horizontal velocity in these 
beds as computed from the short refraction weathering profiles was 
usually close to this value. 

The horizontal velocity in the syenite as determined by plotting 
the times of arrival against the receptor distances for shots in reverse 
directions, and eliminating the apparent velocities due to dip of the 
refracting beds, varied from 10,000 to 25,000 feet/sec. This range in 
velocity was observed in places where syenite was known to exist 
from the core hole data. Due to this large variation in the horizontal 
velocity of the syenite it was impossible to distinguish between the 
syenite and the Paleozoic rocks which have a range of horizontal 
velocities which overlaps the range for syenite. 

Since the survey was made with the idea of obtaining a rapid re- 
connaissance of the area, the control was rather widely spaced except 
at points of unusual interest. It was possible, however, to make a 
map of the Syenite-Paleozoic surface and by using this map in con- 
junction with the map of the top of the Midway derived from the 
core hole data an area could be outlined, outside of which, and lying 
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within the area covered by the survey, it is improbable that deposits 
of commercial bauxite could be found. Within the favorable area it 
is not possible to differentiate between syenite and Paleozoic rocks, 
nor is it possible to predict the presence or absence of bauxite. 

It is regretted that it is not possible to exhibit at this time the 
contoured map of the data, since this map still retains considerable 
commercial value to the company for which the survey was performed. 
It is possible, however, to obtain some idea as to the reliability of the 
data by referring to Table I which gives several comparisons between 
the seismic data and the core hole data. In general, the seismic data 


TABLE I. COMPARISON OF CORE AND SEISMIC DATA 











Hole Formation Core Seismic 
No. by Core Depth Depth 
A-8 Limestone 32 226 
A-9 Limestone 75 273 

C-1-A Syenite 310 321 
C-1 Syenite 225 228 
C-2 Syenite a6 103 
C-3 Syenite 172 173 
C-4 Limestone 272 340 
C-13 Syenite 227 258 
D-1 Syenite 148 178 
F-1 Syenite 145 149 
F-2 Syenite 137 167 





give greater depths than the core hole data, varying from 1o feet 
at core hole C-1-A to 30 feet at core hole F-2. Part of this discrepancy 
is due to the fact that the true first arrival of energy was not read on 
the records, but instead the first trough of the refraction was read and 
used in the computations. This was done since it was thought that 
this trough could be read with more accuracy than the true first 
arrival of the energy. If all the seismic depths are raised 20 feet to 
correct for this difference it is seen that an accuracy of 20 feet is 
indicated. It is probable that some of this variation is due to the fact 
that several of the old core hole locations could not be identified in 
the field and had to be approximated. Coring subsequent to the com- 
pletion of this survey seems to indicate a slightly greater accuracy. 

The authors wish to express their appreciation to Mr. O. C. 
Schmedeman of the Reynolds Mining Corporation for permission to 
publish this paper, and to Dr. George C. Branner, Arkansas State 
Geologist. 





NOTE ON THE TRANSMISSION OF RADIO WAVES 
THROUGH THE EARTH* 


DANIEL SILVERMAN anp DAVID SHEFFETt{ 


ABSTRACT 


_ Experiments are described and data reported on the attenuation suffered by radio 
signals in passing through shallow layers of the earth. Comparison is made between 
the measured attenuations and values computed from published formulae. 


Various investigators have proposed, from time to time, methods 
of geophysical prospecting wnich involve as part of their mechanism 
of operation the transmission of radio waves through the earth: Each 
time such methods are proposed, they revive again the question of the 
possibility of transmission of radio waves through the earth. Although 
a number of studies have been reported upon in the literature,’ the 
question has not been entirely closed. The reason for this, of course, 
lies in the fact that the properties of the earth vary over such a wide 
range and the effect on radio transmission of resistivity, permeability, 
dielectric constant, anomalous dispersion, etc., of the earth over the 
large range in frequencies possible, is not known to a sufficient degree 
to justify a direct statement in all cases. 

Most authorities in the literature provide data which indicate 
that no great penetration of radio waves in the earth is possible. How- 
ever, transmission is known to be a function of the properties of the 
earth and of the frequency and, in general, improves with increase in 
earth resistivity and with decrease in frequency. For all but the lowest 
frequencies the attenuation is considered to be high enough to assign 
a very low probability to the transmission of radio waves to any useful 
depth. 


* Read by title at the 1942 Spring meeting, Denver, Colorado. 

} Joint Exploration Laboratory of the Stanolind Oil and Gas Company and 
Western Geophysical Company. 

t Naval Ordnance Laboratory, Washington, D. C. 

1 Morecroft. Principles of Radio Communication. Wiley & Sons, N. Y. 2nd Ed. 
1927. p. 844. 

Eve & Keys. Applied Geophysics. Cambridge University Press. 1929. pp. 136, 221. 

Hund. Phenomena in High Frequency Systems. McGraw-Hill. 1936. p. 335. 

Nettleton. Geophysical Prospecting for Oil. McGraw-Hill. 1940. p. 380. 

Heiland. Geophysical Exploration. Prentice Hall. 1940. p. 811. 

Feldman. Optical Behavior of the Ground for Short Wavelengths. Proc. I.R.E., Vol. 


21, 1933. P- 293. 
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In contrast with the pessimistic nature of these data are experi- 
ments such as those reported? in which radio signals were received 
in caverns and mines after transmission through considerable thick- 
nesses of rock. Unfortunately, complete data showing the resistivities, 
frequencies and field strength measurements on and below the sur- 
face are not given, so that the actual attenuation of the signals can 
only be surmised. The apparent exceptions provided by these ex- 
amples have kept the question open and have undoubtedly inspired 
similar though unreported measurements of the transmission of radio 
waves through the ground to mines and caves. 
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Fic. 1. Plan and elevation of test location, showing azimuths 
and distances of the stations observed. 


An opportunity was recent! presented to make some measure- 
ments of the attenuation of radio waves in the earth which are here 
reported in the hope that these data may throw some additional light 
on the subject. 

The transmission measurements were made in an active coal mine 
near Tulsa, Oklahoma, at depths of 65 and 75 feet below the surface 
and at 50 and 1100 feet horizontally from a vertical entrance shaft. 


* Eve, Keys, Lee. Proc. I.R.E., Vol. 17, No. 11, 1929, pp. 2072-74. 
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The location is shown in plan and elevation in Fig. 1. In the absence 
of a suitable transmitting source, and also to avoid the difficulties 
of holding field intensity constant very close to a radiating antenna, 
it was decided to use the transmission from regular broadcasting sta- 
tions. Field strength measurements were made with an RCA type 
75B field intensity meter. 

To obtain significant results it was necessary to select broadcast- 
ing stations whose field intensities were constant in this region. This 
ruled out measurements on distant stations. However, enough 
variation in distance of station was provided to show that this factor 
was not of too great importance. It was also decided to measure sta- 
tions distributed in frequency throughout the broadcast band in order 
to note the effect of frequency on attenuation. Selection of stations 
also involved choice of those ranged in different azimuth directions so 
that assurance could be had that the field strength measured was due 
to the wave coming directly from the transmitter. The stations se- 
lected to meet these requirements were KGGF, 690 kc., at a distance 
of 75 miles North 10° East, KVOO, 1170 kc., at a distance of 7 miles 
South 85° East, and KTUL, 1430 kc., at a distance of 5 miles North 
30° West. The directions of these stations with respect to the measur- 
ing point are shown in Fig. 1. Measurements on the surface showed 
that the fields of the three stations selected were uniform in this area. 

In the horizontal shaft of the mine where the measurements were 
made a 220 volt A.C. power line extended for a distance of 1,000 feet. 
There was also a 2 inch pipeline along the side of the shaft which was 
well grounded and extended more than 1100 feet from the vertical 
entrance shaft. Measurements were made in this horizontal shaft at 
two stations. Station A was 50 feet, and Station B, 1100 feet from the 
vertical entrance. Station B was therefore 100 feet beyond the end 
of the power line but was only 5 feet to the side of the pipeline. Station 
A was 65 feet and Station B 75 feet below the surface. 

The data shown in Table 1 indicate the attenuation of the signals 
measured at stations A and B while Table 2 shows the directions from 
which the signals were received. At station A the signal strength was 
reduced by about 30:1 for each of the transmitters without any ap- 
parent relationship between attenuation and frequency. At this loca- 
tion the plane of the loop antenna for maximum signal was always 
directed North-South regardless of the station to which the receiver 
was tuned. In other words, the loop had to be pointed towards the 
vertical mine shaft for maximum signal strength. 
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TABLE 1. FIELD STRENGTH OF BROADCAST STATIONS 











Field Strength 
Broadcast Freq. OnSurface Station A Station B 
Station KC iD Sa fo Ss/fa fe/fo 
KGGF 690 2.1 .063 .O12 33-4 175 
KVOO II70 123 4533 .148 28.4 833 
KTUL 1430 60.6 2.02 .040 30.0 1500 





Since the attenuation at station A was considerably less than at 
station B, even though both points were roughly the same distance 
below the surface, and since the direction from which the signals ap- 
proached the receiving instrument were the same for each signal, and 
exactly the direction of the shaft, the conclusion is indicated that all, 
or at least a major portion of the energy measured at station A ar- 
rived by way of the vertical and horizontal shafts and power line, 
rather than directly through the ground. The results obtained at this 
station can therefore be disregarded. 

At station B the direction of the loop for maximum signal strength 
was in each case very close to the true direction of the broadcast sta- 
tion tuned in. The discrepancies in the cases of KGGF (10°) and 
KTUL (4°) fall within the limits of error due to the broad indication 
of the loop. In the case of KVOO the 18° difference is not easily ex- 
plained on this basis, and indicates that the resultant signal received 
at station B is made up of two parts, one coming directly from the 
broadcast station through the earth, and another part arriving at B 
by way of the vertical and horizontal shafts. A simple vector compu- 
tation indicates the direct component to be approximately 92% of 
the total. It can be said, therefore, that at least the major portion 
of the energy of the broadcast signals had to penetrate the earth in 
order to reach the loop antenna of the measuring set. 

The attenuation of signals at this location was much greater than 


TABLE 2. DIRECTIONS OF PLANE OF LOOP 
FOR MAXIMUM SIGNAL 











ae mg On Surface Station A Station B 
KGGF N 10° E N N 
KVOO S 85° E N S 67°E 
KTUL N 30° W N — 





Bearings are very broad, being good to about +7 degrees. 
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Fic. 2. Attenuation of radio frequency waves after passage through 75 feet of 
earth. Circles are observed values, crosses are computed on the basis of a resistivity 
of 5500 ohm-cm., triangles are computed on the basis of 12,000 ohm-cm. 
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that near the mine entrance, and also exhibited a distinct variation 
with frequency as shown in Table 1. These data are plotted in the 
upper curve of Fig. 2. 

In order to compare the measured results with the theoretical 
attenuation of radio signals in the earth as given by A. S. Eve,’ it was 
necessary to have the value of resistivity of the earth in the region 
of the mine. A shallow electrical resistivity depth profile, interpreted 
in accordance with the curves of Wetzel and McMurray,‘ gave the 
resistivity of the top layer as 12,000 ohm-cm and the thickness of 
this layer as 105 feet. The resistivity of the intermediate layer was 
indicated to be 1200 ohm-cm. and its thickness 35 feet. 

The relationship given by Eve shows the attenuation factor to be 


- [- an {* (e " ay" 7 ‘ " 


where f=frequency in. cycles per second 
d=depth in centimeters 
c=velocity of light 
u=permeability of the earth 
e=dielectric constant of the earth divided by 2 
o=conductivity in electrostatic units 





Using this formula and the value of resistivity of 12,000 ohm-cm, the 
attenuations for the three frequencies observed are plotted in Fig. 2 
in the lower curve. A very large discrepancy is noted between observed 
and computed values of attenuation. The observed attenuation for the 
case of KVOO, for instance, was approximately ten times as great as 
the computed value and similarly, for the other two stations, the 
attenuations observed were many times those computed on the basis 
of an earth resistivity of 12,000 ohm-cm. The general shape of the 
curve provided by the computed data was similar to that of the ob- 
served data and an attempt was made to choose a different value of 
resistivity such that it would make the attenuation corresponding to 
the frequency of 1430 KC the same in both cases. This condition was 
satisfied by a value of resistivity of 5500 ohm-cm. Using this value, 
the attenuations corresponding to the other two frequencies were com- 


3 A. S. Eve. Technical Publication No. 316, Class L, Geophysical Prospecting, No. 
19. A.I.M.E. Presented at New York meeting February, 1930. 
4 Wetzel & McMurray. Three Layer Resistivity Problem. Geopxysics, Vol. 2, 


No. 4, (Oct. 1937), pp. 321-341. 








DANIEL SILVERMAN AND DAVID SHEFFET 


s000 és 





1000 ra 





400 








100 1 


7 ae 











hl 


























if 





° goo 800 1200 1600 2000 
Ke 
Fic. 3. Attenuation of radio frequency waves after pasasge through 75 feet of 
earth, extended frequency range. Circles are observed values, solid curve is computed 
on the basis of a resistivity of 5500 ohm-cm. 


RS AO IN Lies Al aint soir a IRS SE I a * 






by SARA ete ts DEIN RIOR ICES SO cot 








TRANSMISSION OF RADIO WAVES THROUGH THE EARTH 413 


puted and are shown plotted in Fig. 2 superimposed on the observed 
curve. Using this value of resistivity, all three points follow closely 
the curve obtained by actual measurement of field strengths in the 
mine. 

A rough substantiation of this value of resistivity is found in the 
data given by the Federal Communications Commission on ground 
conductivity in northeastern Oklahoma (F.C.C. Ground Conductivity 
Map of June, 1938). These data are derived from broadcast field 
measurements on the surface of the ground. The ground conductivity 
in the region of Tulsa is given as 15 X10-“ e.m.u., which corresponds 
to a value of resistivity of 6600 ohm-cm. This is comparable with the 
value of 5500 found necessary to make a direct correspondence be- 
tween observed and computed attenuations. 

These experiments are admittedly incomplete since they simply 
give an over-all figure of loss and do not indicate what portion, if 
any, of the total loss occurs in the transition of the wave across the 
boundary plane of air and earth at the surface, and the similar plane 
in the tunnel, and what portion is due to the absorption in the 75 feet 
of earth. Assuming that the measured attenuation is truly representa- 
tive of the attenuation per unit of distance in the earth, a further step 
has been taken using Eve’s formula to extend the attenuation curve 
to higher and lower frequencies. This is shown in Fig. 3. It is of interest 
to note that the theoretical attenuation at this location at 50 KC is 
only 4:1 for a depth of 75 feet as compared to 1500:1 at 1400 KC. 
However, while the attenuation is very considerably reduced for low 
frequencies, the practical difficulty of radiating energy is almost in- 
surmountable, requiring for a quarter wave length antenna at 50 KC 
a wire of approximately 1500 meters in length. 

This report is admittedly incomplete, in many respects, and par- 
ticularly as to the effect on attenuation of the proximity of the trans- 
mitter, and the directivity of the antenna system. This was, of course, 
impossible to determine under the experimental conditions. 

Appreciation is expressed to Mr. L. Watt Stinson and to the 
management of station KVOO for the loan of the field strength meter 
used in making these observations, and to the officers of the Stanolind 
Oil and Gas and Western Geophysical Companies for permission to 
publish this note, 











THE IMPEDANCE OF A GROUNDED WIRE* 
ALFRED WOLFt 


ABSTRACT 


The impedance of an insulated wire stretched along the surface of the earth, 
regarded as a homogeneous conductor, is a function of frequency and of the con- 
ducitivity of the earth. Formulas are given for the inductance and the resistance of 
such a wire which are applicable under conditions met with in geophysical prospecting. 


In connection with some problems of electrical prospecting, it is 
sometimes required to calculate the inductance of a long insulated 
wire stretched along the surface of the earth and grounded at both 
ends. The earth, which is regarded as a semi-infinite homogeneous 
conductor, forms the return part of the electric circuit and causes the 
impedance to vary with frequency. This is the case because the cur- 
rent penetration into a conductor is a function of frequency, and 
determines both the inductance and the a.c. component of the re- 
sistance of the grounded wire. The d.c. component of resistance is a 
function of the size of the wire and of the electrodes used for ground- 
ing, and will not be considered further. 

In the following, it will be assumed that the length of the wire is 
great compared to its diameter and to the dimensions of the ground- 
ing electrodes. Further, the frequency is taken to be low enough, so 
that the skin effect in the wire itself is negligible. These conditions 
are valid for a wide range of circumstances met with in the practical 
application to gephysical prospecting. 

Carson! has solved a problem similar to the above on the assump- 
tion that the length of the wire is great compared to the wave length; 
his results can be employed in calculations dealing with long tele- 
phone lines. The fundamental assumption made in the following is 
that the current along the grounded wire is constant. The calcula- 
tion is based on formulas given by Foster? for the mutual impedance 
of grounded wires. The earth is regarded as a homegeneous con- 
ductor of negligible magnetic susceptibility. 

Let a unit current of frequency f be passed through the wire. The 
integral of the electric field along a straight line parallel to the wire, 


* Read by title at the Denver meeting, April, 1942. 
+ Geophysical Research Corporation, Tulsa, Oklahoma. 

1 J. R. Carson: Bell System Technical Journal 5, 1926, p. 530. 
2 R. M. Foster: Bell System Technical Journal 10, 1931, p. 409. 
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at a distance x from it, and of the same length as the wire, is, accord- 
ing to Foster, given by 
M(x) = F(x)/2md (1) 


where 





F(a) = 3 do i, ac . — (2) 


with 





r= vV/x2+ (u — 2)? 
y = Vanjor 





w = 2nf 

javan 

2a = | = length of wire 

dX = conductivity of the earth in e.m.u.* 


Since the distribution of current in the wire is uniform, the a.c. 
component of impedance will be obtained approximately by aver- 
aging M(x) over the cross-section of the wire. This step is necessary, 
as in all cases of self-inductance calculation, because M(x) becomes 
logarithmically infinite as « approaches zero. The total impedance of 


of the circuit is 
Z = R+F(x)/20r (3) 


where R is the d.c. resistance in e.m.u. The resistance of the wire itself 
is included in R. 
Since the radius of the wire is small compared to its length, F(x) 
has to be determined only for small values of x. 
If the variables 
I 
— (u — 2) 


et 


(4) 


I 
n=—(ut9) 
V2 


* All quantities are expressed in absolute electro-magnetic units. 1 Ohm=1o0® 
e.m.u., 1 Henry=10% e.m.u. 













416 ALFRED WOLF 


are introduced, the integrand in (2) becomes independent of 7. The 
integration with respect to 7 introduces the factor. 


2 za| 1 - Le | 
V2a 





and after substituting 














V2t=t 
(a)/a= [at ae [1H] 
Ez f 1—(1+7V 22+?) exp [—yV/ 2? +? | (5) 
——| tdi | 
l 0 (x? £?)3/2 





From (s) it follows easily 






l ile 
Mp as (1—exp [—yvV/x?+/?]) 


—2 ff exp [--WF FF lu (6) 
vd 0 







I I 721 ]2) 
a (x et (1—exp [—vV2?+?]). 






Since only terms of the lowest order in x are to be preserved this can 
be written . 
I — ev! ~ cameeani 
F(a)/2t =~ —S— — 7 [ dwexp [- WF FH] (0) 
x? x? J 9 
v r—e ve ye vt 











2)? 


al 











Now, since 


fi aver [— yu? + 27] -[-f- «Ki(yx) -f 





and since 










yx? yx I : 7; 
4+Ki(yx) = — log — + — + terms of higher order in % 
2 2 7 
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it follows that 


7 2 ° du 
F(x)/2l = ~ | tog = — f — on | 
2 yx 1 4Uu 





(8) 
I 
+ — [1 — e+ ye”! — ay]. 
al? 
Averaging log 1/x over the cross-section of the wire gives 
I I I 
log — = log — + — (9) 
x b 4 
where 3 is the radius of the wire. 
Combining (3), (8), (9) it follows that 
2 I ° du 
Z=R + ja:at| log + — -f —¢ 
yb 4 yl U 
(10) 


I 

+ yp (rx —eve+ lll _ oat]. 
The evaluation of the complex logarithmic integral in (10) would 

involve considerable labor since tables are not available. Fortunately, 

the formula can be simplified both for low and for high frequencies. 
When y/<1 


° du 
f —e-™ = — log 1.7817l.* 
1 U 


It follows that, for low frequencies, (y/1), when the current pene- 
tration is greater than the length of the wire 





; 3.5621 
Z=R + jell log a 1.25 (11) 


and the low frequency inductance is 


l 
ll tog " + 02 | , 


When the current penetration is smaller than the length of the 
wire, so that yi>1, but still yb<z, the exponential terms in (10) 
can be neglected and ; 


* Cf., Jahnke-Emde: Tables of Functions, 2nd ed., p. 79. 
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Z=R+ Fae it I +2] 
= a jw2 ee | (12 


The expression 





I I 
il log a =| (13) 
for the high frequency inductance of the circuit decreases slowly as 
the frequency increases. 

A certain amount of resistance is introduced into the circuit be- 
cause of the concentration of current near the surface of the ground. 
This a.c. component of resistance is given by the second term on the 
right hand side of (12) 

T 
— wl. (14) 
4 

As an illustration, the following problem may be considered. A No. 
18 copper wire conductor (b=0.051 cm) a mile long (J= 1.609 X 10° cm) 
is stretched along the surface of the earth which is assumed to have 
a resistivity of 1000 ohm-cm (A=107"* e.m.u.). It is desired to calcu- 
late the resistance which is needed in the circuit in order that its 
time constant be 0.0001 sec. The inductance will be calculated for a 
frequency of 2500 cps, since that makes the duration of a quarter 
cycle equal to the time constant. The value of | y/| at this frequency is 
71.4, and, hence, the inductance is obtained from (13), and the high 
frequency ground resistance from (14). An easy calculation shows the 
inductance to be 0.00374 henrys and the a.c. component of ground 
resistance, neglecting the resistance at the grounding electrodes, is 
2.0 ohms. A time constant of 0.0001 sec. requires a total circuit resist- 
ance of 37.4 ohms and, consequently, the value of the d.c. resistance 
R should be 35.4 ohms. This happens to be very nearly the value of 
the resistance of the copper conductor itself which is 34.4 ohms at 
<<. 











DISCUSSIONS AND COMMUNICATIONS 


To the Editor: 


In his paper in the April and July 1942 issues of GEopHysics, Dr. Sharpe gives a 
very enlightening analytic solution to the problem of the wave motion generated in an 
elastic medium by application of a pressure of arbitrary form to the interior surface 
of a spherical cavity. It is very encouraging to find that deductions regarding the rela- 
tion of the amplitude and frequency of the generated wave to the physical character- 
istics of the medium agree so well with empirical operation. A thorough understanding 
of the explosion process and its accompanying wave phenomena is a prerequisite to 
the final refinement of the seismograph as a tool of Exploration Geophysics. It is to 
be hoped that Dr. Sharpe and his colleagues, as well as others, will continue the com- 
bined analytic and experimental attack. 

In connection with this discussion of Dr. Sharpe’s paper, I wish to call attention 
to A.I.M.M.E. Technical Paper 1268, in the Petroleum Division Transactions for 
1941, by Evinger and Muskat, on “Some Theoretical Considerations on the Problem 
of Well Shooting,” that is, on the problem of shooting oil wells for the purpose of 
increasing production. To a considerable extent these papers are complementary. Dr. 
Sharpe is principally concerned with the wave motion outside the sphere or zone of 
fracture, whereas, Dr. Muskat and Mr. Evinger are concerned with the stresses in the 
medium and with the interior of the sphere or zone of fracture. Geophysicists, by virtue 
of their training and instrumental facilities for measuring displacement, velocity, 
acceleration, and stress in a medium, and density, seismic wave velocity and elastic 
parameters of rocks, are eminently fitted to improve our knowledge of this important 
production problem. 

GeEorcE B. LAMB 
Houston, Texas, 
August 26, 1942. 


REVIEWS 


A Magnetic Survey of Central South Dakota. Geophysics by W. H. Jordan. Geology by 
E. P. Rothrock. South Dakota Geological Survey Report of Investigations No. 37 
(Vermillion, Nov., 1940). 35 pp., 8 illus., including magnetic map of area. 


Magnetometer Surveys During 1941, by Edward L. Tullis. South Dakota Geological 
Survey. Report of Investigations No. 42 (Vermillion, March, 1942). 40 pp., 13 
illus., including magnetic map of area. 


These are the two latest reports on the program of magnetic observation which has 
been undertaken by.the South Dakota State Geological Survey. According to R. I. 
No. 37 the plan is to make a magnetic survey of the entire state. The purpose of this 
program according to R. I. No. 37 is “first to obtain information for the people of the 
state, and second to encourage further activity on the part of the oil companies.” To 
accomplish this purpose they are making a regional magnetic (vertical component) 
map of the state. Magnetic stations are located approximately six miles apart. Most 
of the work to date covers a large area located in the central part of the state and 
extending from the southern to the northern boundary. R. I. No. 37 presents the verti- 
cal magnetic data at 118 stations, covering approximately 10,000 square miles. R. I. 
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No. 42 gives the same data for 237 additional stations, covering an area of approxi- 
mately 6850 square miles. 

Both reports discuss very briefly the theory and field technique of magnetic sur- 
veying as well as the construction of the magnetometer. Both reports also discuss 
briefly the theory of interpretation of magnetic data, and offer some interpretation 
of the data obtained by these two surveys. 

All station values have been corrected for difference in latitude and longitude. The 
exact corrections used are given. A map showing the amount of the error of misclosure 
is included in each report. The map of misclosure, the discussion of errors, and the dis- 
cussion of the field methods which were used give one a good idea of the accuracy of 
the work. 

A contoured magnetic map is included in each report. The contour interval used 
is 100 gammas. The station locations are shown on the map, but not the values at the 
individual stations. The individual station locations and values are tabulated so that 
the data may be transferred to any other map for additional study. Any geologist 
or geophysicist interested in South Dakota will welcome the activity of the State 
Geological Survey in preparing these reports. 

GarvIN L. TAYLOR 


Seismology, by Perry Byerly. 256 pp. Prentice-Hall, New York, 1942. 


The book is divided into two parts, the first of which appears to be lecture notes 
for a popular course in seismology. The various theories and speculations as to 
the immediate and underlying causes of earthquakes are set forth and the effects of 
earthquakes are described with numerous photographic illustrations. Brief descriptions 
of great earthquakes of history are given, the earliest one described being the St. 
Lawrence Valley earthquake of 1663. In a chapter on the distribution of earthquakes, 
the author assigns a seismicity of 10 to the San Francisco Bay Area and 3.6 to the Los 
Angeles Area and one wonders how long these figures will remain unchallenged. 

Part II of the book covers instrumental and wave theory, instrumental observa- 
tions, and interpretations. Seismograph theory is developed using operational methods. 
Elastic wave theory is set forth in almost outline form. Seismic prospectors should find 
the chapters on “Paths of Waves and Travel Times” and “Location of Epicenters” of 
interest. In reading them, they should bear in mind that the seismologist ordinarily 
works without a time break, with no exact knowledge of the location of the “shot 
point,” no weathering information, no well velocity surveys, and no correlation 
spreads. However, the simplicity of having to deal with only one recognized sub- 
surface reflecting interface, may seem adequate compensation. 

The author has achieved his desire for brevity and provides references at the end 


of each chapter for those interested in detail. 
Kari Dyk 























PATENTS* 


ELECTRICAL PROSPECTING 


U.S. No. 2,284,990. M. Schlumberger. Iss. 6/2/42. App. 7/5/39. 


Method and A pparatus for Electrical Prospecting, A method of electrical prospecting 
making use of telluric earth currents in which the earth potentials or their derivatives 
are recorded by frequency selective devices. 


U.S. No. 2,288,310. T. Zuschlag. Iss. 6/30/42. App. 4/23/38. Assign. Lundberg Ex- 
ploration S. A. 


Apparatus for Geoelectric and Seismic Investigations. A method of comparing steady 
state a-c voltages obtained either in a-c electrical prospecting or continuous wave 
seismic prospecting by using two amplifiers whose outputs are connected in a bridge 
circuit and balanced. 


GEOCHEMICAL PROSPECTING 
U.S. No. 2,284,147. H. N. Herrick. Iss. 5/26/42. App. 1/28/39. Assign. Standard Oil 
Co. 


Method and A pparatus for Soil Gas Analysis. A method of collecting soil gas samples 
by sealing a cavity in the soil, permitting atmospheric breathing, and absorbing hydro- 
carbon gas from the air passing in and out of the cavity during the breathing process. 


U.S. No. 2,287,059. R. N. Platts, F. B. Slichter, E. L. O’Donnell, and P. T. Bennett. 
Iss. 6/23/42. App. 5/25/40. 
Apparatus for Obtaining Soil Samples. A sample or core taking machine which 
forces a hollow cutter into the earth by hydraulic pressure. 


U.S. No. 2,287,101. L. Horvitz. Iss. 6/23/42. App. 1/8/38. Assign. E. E. Rosaire. 


Means and Method of Analysis. A method of analyzing soil gas for minute traces 
of hydrocarbons in which the gas is circulated repeatedly around a closed loop having 
a refrigerating and a combustion zone. 


SEISMOGRAPH PROSPECTING 


U.S. No. 2,285,610. O. S. Petty. Iss. 6/9/42. App. 8/1/40. 


Method and Apparatus for Seismic Surveying. A seismograph amplifier which dif- 
ferentiates the initial wave form. 


U.S. No. 2,286,106. O. F. Ritzmann. Iss. 6/9/42. App. 2/9/39. Assign. Gulf Research 
& Development Co. 


Apparatus for and Method of Receiving and Recording Vibrations. A seismograph 
amplifier having an automatic characteristic control which is effected by the first 
arrivals but which is subsequently cut off and thereafter acts as a function of time 
only. 


U.S. No. 2,286,170. C. A. Heiland. Iss. 6/9/42. App. 7/26/37. Assign. Heiland Research 
Corp. 


* Abstracts by O. F. Ritzmann, Gulf Research & Development Co. 
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Selector System for Amplifiers. A system of controlling the amplification of seismo- 
graph amplifiers according to any predetermined function of time, control being ac- 
complished by means of a time sequence of shutters interposed between a lamp and a 
photocell which effects the amplification. - 


U.S. No, 2,286,386. D. Silverman. Iss. 6/16/42. App. 12/7/39. Assign. Stanolind Oil 
and Gas Co. 
Seismometer, A suspended magnet type of dynamic geophone having two elements, 
the coil of each being rigidly attached to the magnet of the other and the coils elec- 
trically connected so as to eliminate response to rotational motion. 


U.S. No. 2,286,567. J. O. Parr, Jr. Iss. 6/16/42. App. 9/12/40. Assign. Olive S. Petty. 

Seismic Surveying. A vacuum tube circuit for combining the outputs of a number 
of geophones so that response is blocked unless all geophones are delivering voltage of 
the same sign. 


U.S. No. 2,286,897. P. J. Costa and G. E. Holback. Iss. 6/16/42. App. 12/20/38. 
Vibration Pickup. A dynamic type vibrometer having a suspended magnet and 
stops connected so as to short circuit the coil if the vibrating amplitude becomes ex- 


cessive. 


U.S. No. 2,288,838. .. W. Pike and A. W. Vance. Iss. 7/7/42. App. 2/29/40. Assign. 
Radio Corp. of America. 
Vibration Pickup Device. A displacement type vibrometer in which the vibration 
changes either the capacity or the inductance in a resonant circuit connected to a 
source of constant frequency a-c slightly off resonance. 


WELL LOGGING AND SURVEYING 


U.S. No. 2,284,345. C. H. Schlesman. Iss. 5/26/42. App. 10/21/41. Assign. Socony- 

Vacuum Oil Co. 

Method and Apparatus for Geophysical Prospecting. A gamma ray well logging 
instrument using a number of detectors which are connected in sequence to the record- 
ing mechanism, the sequence rate being related to the rate of movement of the instru- 
ment so as to give improved resolving power. 


U.S. No. 2,285,809. R. F. Davis. Iss. 6/9/42. App. 4/4/40. Assign. Well Surveys, Inc, 

Well Surveying Method and Apparatus. A method of transmitting information ob- 
tained by a well logging device to the surface by means of mechanical vibrations sent 
up through the supporting cable. 


U.S. No. 2,285,840. S. A. Scherbatskoy. Iss. 6/9/42. App. 8/4/39. Assign. Well Surveys, 

Inc. 

Well Survey Method and Apparatus. A gamma ray well logging device using two 
ionization chambers, one of which responds to all radiations and one shielded to respond 
only to hard radiations, and recording the difference in their ionization currents in 
order to eliminate the effect of sources close to the device. 
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U.S. No. 2,288,278. L. G. Howell. Iss. 6/30/42. App. 6/23/41. Assign. Standard Oil 

Development Co. 

Gamma Ray Well Logging. A method of obtaining gamma ray logs from an instru- 
ment having a number of counters but using only a single cable by filtering the pulses 
from each counter before impressing them on the cable and again filtering at the sur- 
face to identify the pulses from the respective counters. 


U.S. No. 2,288,876. J. C. Arnold. Iss. 7/7/42. App. 1/21/41. Assign. Lane-Wells Co. 

Magnetic Logging. A magnetic well logging device in which an inductor coil is 
rotated simultaneously about a horizontal and a vertical axis while traversing the hole 
and recording the instantaneous form of the emf generated. 


U.S. No. 2,288,884. M. C. Bowsky. Iss. 7/7/42. App. 12/11/39. Assign. Lane-Wells 

Co. 

Electrical Logging in Oil-Filled Wells. An a-c electrical well logging method in which 
the electrodes make contact with the walls of the hole and using a constant current 
a-c source, the pick up circuit being inductively coupled to the exciting circuit so as 
to eliminate contact resistances. 


U.S. No. 2,288,973. J. Neufeld and E. H. Cooley. Iss. 7/7/42. App. 12/15/39. Assign. 

Well Surveys, Inc. 

Well Surveying Method and Apparatus. A method of determining dip of under- 
ground strata traversed by a well by observing the asymmetry of the current field 
when one terminal of a battery is connected to an electrode in the well and the other 
terminal is connected through resistors to two surface electrodes diametrically opposite 
and equidistant from the well. 


U.S. No. 2,289,926. J. Neufeld. Iss. 7/14/42. App. 3/18/40. Assign. Well Surveys, Inc. 

Well Survey Met'od and Apparatus. A method of radioactivity well logging in 
which a record is made of the natural radioactivity and another record of the com- 
bined natural plus induced radioactivity, and by subtraction arriving at a record of 
the induced radioactivity alone. 


U.S. No. 2,290,075. M. Schlumberger. Iss. 7/14/42. App. 4/3/39. Assign. Schlumberger 

Well Surveying Corp. 

Thermal Process and Device for Surveying the Beds Traversed by Drill Holes. A 
method of obtaining a thermal property log of a well by simultaneously heating an 
extended portion of the well fluid by a current carrying resistance wire and making a 
temperature survey before and after the heating. 


MISCELLANEOUS 


U.S. No. 2,284,654. B. M. Harrison. Iss. 6/2/42. App. 10/20/38. Assign. Submarine 

Signal Co. 

Distance and Depth Measuring Device. An echo-sounding device in which correc- 
tions for water temperature and salinity are automatically made by the method of 
having the depth indicator driven through a differential gearing by two adjustable 
speed electric motors. 
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U.S. No. 2,286,673. L. A. Douglas. Iss. 6/16/42. App. 6/10/41. 


Means for Extracting the Core Content of Subterranean Strata. A side wall sampling 
device having a mechanically operated plunger which is forced into the side wall and 
a vacuum chamber which takes up gas or liquid from the side wall cores. 


U.S. No. 2,287,819. P. B. Nichols. Iss. 6/30/42. App. 5/21/38. Assign. One-half to 
Cities Service Oil Co. 
Device for Recording Drilling Operations. A rate-of-drilling recorder which records 
time, depth, speed of penetration, vertical chatter of drill pipe and other drilling 
operations. 


U.S. No. 2,288,210. M. Schlumberger. Iss. 6/30/42. App. 5/3/39. 


Core Taking Projectile. An easily withdrawable side wall coring projectile having 
a flange or collar which makes a hole larger than the main part of the projectile and 
which is loosened from the projectile during impact. 


U.S. No. 2,288,942. R. L. Doan. Iss. 7/7/42. App. 6/11/40. Assign. Phillips Petrol- 
eum Co. 
Fluoroscope. A portable apparatus for detecting crude oil in drilling mud in which 
the mud is continuously circulated under a source giving only ultraviolet light and the 
fluorescence is observed through an eyepiece shielded from direct rays of the source. 


U.S. No. 2,289,961. J. E. Hancock. Iss. 7/14/42. App. 5/6/39. Assign. General Electric 
Co. 
Vibration Generator. A dynamic type vibrator for vibration testing of materials 
including a smaller dynamic type pickup device for determining the actual motion 
developed. 


U.S. No. 2,289,962. J. E. Hancock. Iss. 7/14/42. App. 5/6/39 and 11/6/41. Assign. 
General Electric Co. 


Vibraton Device. A dynamic type vibrator having a locking device to prevent 
straining the moving system when specimens are being mounted on it. 
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D. C. SKEELS received a B.A. degree from 
the University of Montana in 1930, with a 
major in physics and mathematics, after which 
he spent three years as a Rhodes Scholar at 
Oxford, where he studied geology. He then 
studied for the Ph.D. degree in geology at 
Princeton, receiving his degree in 1936. Dr. 
Skeels has been engaged in geophysical work 
since 1936, in Tulsa for the Carter Oil Com- 
pany, and in New York City for the Standard 
Oil Company of New Jersey. He is a member 
of Sigma Xi, the American Association of Pe- 
troleum Geologists, the American Geophysical 
Union, and the Society of Exploration Geo- 
physicists. 





D. C. SKEELS 





Brrp G. Swan received his B.S. degree in 
Geological Engineering from the University of 
Oklahoma in 1932, after having spent the two 
previous years with Continental Oil Company 
in torsion balance and seismograph field work. 
Since 1932 he has been associated with Con- 
tinental Oil Company as seismograph com- 
puter, party chief, and since 1937 in inter- 
pretative and administrative work in the 
Ponca City office. Mr. Swan is a member of 
the Society of Exploration Geophysicists. 
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Mitton B. Dosrin was graduated from 
the Massachusetts Institute of Technology in 
1936 and subsequently attended Columtia 
University and the University of Pittsburgh, 
receiving his M.S. at the latter institution in 
1941. In 1937 he became connected with the 
Geophysics Division of the Gulf Research and 
Development Company for whom he was en- 
gaged in seismic interpretation both at the 
company’s Pittsburgh headquarters and in 
the field. In June, 1942, he joined the staff of 
the Naval Ordnance Laboratory, Washington, 
D. C., as a contract physicist. He is a member 
of the American Physical Society, the Ameri- 
can Geophysical Union and the Society of 
Exploration Geophysicists. 
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LorENz SHOcK attended the University 
of Oklahoma, where he received an A.B. de- 
gree in 1937. After a summer with the Wes‘ern 
Geophysical Company he returned to the Uni- 
versity as a research scholar, receiving his 
M.S. degree in 1938. Mr. Shock has been em- 
ployed by the National Geophysical Company 
since 1938, with the exception of one year spent 
at the University of Oklahoma as a graduate 
student and part time instructor in physics. He 
is a member of Phi Beta Kappa and an associ- 
ate member of Sigma Xi. 
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E. D. Atcock received an A.B. degree in 
mathematics from U.C.L.A. in 1930, after 
which he worked for the Bell Telephone Labo- 
ratories for two years as a research engineer. 
He then returned to take his Ph.D. degree in 
physics at the California Institute of Tech- 
nology, in 1935. After receiving his Ph.D., Dr. 
Alcock accepted employment in the geo- 
physical industry, first with the Shell Pe- 
troleum Corporation, for whom he worked 
three years as a seismologist, and then with 
National Geophysical Company, with whom 
he is now employed. 
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DANIEL SILVERMAN received Lis B.S. de- 
gree in 1927 from the University of California. 
He attended the Massachusetts Institute of 
Technology from 1928-30 as the Heller Travel- 
ling Fellow from the University of California. 
There he received his M.S. degree in 1929 and 
Sc.D. in 1930. He was employed as research 
engineer by the Westinghouse Research Lab- 
oratory from 1930 to 1934, and by the Arma 
Corporation of Brooklyn, N. Y., from 1934 to 
1936, where ke worked on the development of 
navigation and fire control instruments for the 
U.S. Navy. He spent the year 1936-37 as Asst. 
Professor of Mechanical Engineering at the 
University of California, and then became as- 
sociated with the Joint Exploration Labo~a‘ory 
of the Stanolind Oil and Gas Company and Dante, Gianni 
Western Geophysical Company in 1937 as 
physicist in charge of design and construction, and since 1940, as laboratory super- 
visor in charge of all research and development activities. 

Dr. Silverman is a member of the American Association of Petroleum Geologists, 
the American Institute of Mining and Metallurgical Engineers, the Society of Explora- 
tion Geophysicists, the Seismological Society of America, the American Geophysical 
Union, the Acoustical Society of America, and the Institute of Radio Engineers. 








Davin SHEFFET attended the California 
Institute of Technology, where he received a 
B.S. degree in electrical engineering in 1930 
and an M.S. degree in 1931, and served as a 
laboratory instructor and took additional 
graduate work in 1932. Mr. Sheffet did elec- 
tronic development work in 1933 for the West- 
ern Electrcnics Company, well-logg’ng re- 
search in 1934 for the Elliott Core Drilling 
Company, and then joined the staff of the 
Joint Exploration Laboratory of the Stanolind 
Oil and Gas Company and the Western Geo- 
physical Company, where he has specialized 
in instrumental developments. Since 1941 
Mr. Sheffet has been on leave from his firm as 
A 4 | a contract employee of the Naval Ordnance 
Laboratory. He is a member of the Society of 
DAVID SHEFFET Exploration Geopkysicists. 
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ALFRED WOLF took his B.S. and M.S. in 
physics at M.I.T. in 1922 and 1923, and after 
four years of work in geophysics at the Mellon 
Institute on behalf of the Gulf Oil Corporation 
and the Geophysical Research Corporation, 
went to Leipzig where he received his Ph.D. in 
1931. Dr. Wolf has been associated with the 
Geophysical Research Corporation since 1935, 
engaging principally in theoretical investiga- 
tions of exploration problems. He is a member 
of the American Association for the Advance- 
ment of Science, and a Fellow of the American 
Physical Society. 
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Otro F. RitzMANN, who prepares the 
patent abstracts which are a regular feature 
of Gropnysics, received his B.S. degree in 
electrochemical engineering from Pennsylvania 
State College in 1922, and his M.S. degree in 
physics in 1924. He continued his graduate 
work in physics at the California Institute of 
Technology. Mr. Ritzmann taught college 
physics at the Pennsylvania State College, 
California Institute of Technology, Michigan 
College of Mines, and Lehigh University and 
in 1929 joined the staff of the Gulf Research 
and Development Co. He is a member of Phi 
Kappa Phi, Sigma Xi, American Association 
for the Advancement of Science, American 
Physical Society, American Geophysical Union, 
Society of Exploration Geophysicists, and the 
Physical Society of Pittsburgh. Otto F. RitzMANN 























ARNOLD J. F. SrEcERT. A biography and photograph will be found in Gropuy- 
sics, Vol. VII, No. 1, (January 1942). 
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ANNOUNCEMENTS 
STUDENT COMPETITION 


At a meeting of the Executive Committee held in Houston on August 26, 1942, 
it was decided that steps be taken to encourage the submittal of papers from graduate 
and undergraduate members of universities. To achieve this it was decided to make 
use of the authority granted to the Executive Committee by the amendment to the 
by-laws adopted at the annual meeting in Denver, under which suitable prizes may 
be awarded for the best papers. A notice announcing these prizes will be circulated to 
all universities where courses in geophysics are customarily given. The proposed form 
of this notice is as follows: 

“The Society of Exploration Geophysicists desires to stimulate study of geophysi- 
cal problems by encouraging the publication by undergraduate and graduate students 
of papers dealing with geophysics. Accordingly, the Society is willing to consider 
such papers for publication in its Journal Gropuysics. In particular it is hoped that 
theses dealing with geophysical problems which have been accepted from candidates 
for degrees will be submitted. 

The Editor, and such committee as he may appoint for the purpose, will referee 
all papers submitted and in the event an outstanding paper is received and published 
which in the opinion of the Editor inerits special recognition, the author may, at the 
discretion of the Executive Committee, be awarded free membership in the Society 
for a period of three years in the class of membership for which he is eligible (i.e., 
Student, Associate or Active). 

Manuscripts should be submitted to the Editor, Joseph A. Sharpe, P. O. Box 591, 
Tulsa, Oklahoma. These should be an original typewritten copy, either double or 
triple spaced, with wide margins. Manuscripts and illustrations of papers not accepted 
for publication will be returned to the authors; manuscripts and illustrations of pub- 
lished papers will be destroyed immediately after publication unless the author requests 
their return.” 


FEDERAL EXPLOSIVES ACT 


Section 18 of the Federai Explosives Act of December 26, 1941, authorizes the 
Director of the Bureau of Mines to issue rules and regulations to effectuate the purposes 
of this Act. The Bureau has issued regulations and appointed explosives investigators 
to supervise the carrying out of such regulations. The supervising engineer for District 
G, embracing most of the territory in the southwestern states where seismograph oper- 
ations are being carried out, is Mr. G. M. Kintz whose office is 544 Federal Building, 
Dallas, Texas. During the Annual Meeting of the Society at Denver, a conference 
with Mr. Kintz was arranged for the purpose of discussing regulations to be promul- 
gated by the Bureau of Mines. A committee of the Society comprising Messrs. G. W. 
Carr, T. I. Harkins, F. F. Reynolds and F. Goldstone was appointed by the President 
to make recommendations to the Bureau of Mines. This committee drafted such recom- 
mendations which were in due course forwarded to Washington. A further meeting was 
held in Dallas at which Mr. Kintz, his inspectors, and representatives of geophysical 
companies were present. At this time further discussion was had concerning the speci- 
fications for field storage magazines and other related matters. 
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Growing out of these various conferences and others held by the Bureau of Mines 
with different branches of industry using explosives, a brochure was issued dated 
July 15, 1942, under the title “Instructions for Storing, Handling and Transporting 
Explosives.” This can be obtained from the Department of Interior, Bureau of Mines, 
Explosives Control Division, at Washington, D. C. These instructions do not lay down 
inflexible rules on all matters, much being left to the judgment of explosives investi- 
gators and other inspectors of the Bureau of Mines in the field; accordingly, the fol- 
lowing letter from Mr. Kintz to Mr. Goldstone under date of September 14, 1942, will 
be of interest. 


Mr. F. Goldstone, President 

Society of Exploration Geophysicists 
P. O. Box 2099 

Houston, Texas 


Dear Mr. Goldstone: 


After the meeting here in Dallas with the representatives of the various geo- 
physical companies, their suggestions were submitted to our Washington office, and 
I am glad to say that the Washington office agreed with what was-decided on at that 
meeting, namely: 


1. Present geophysical magazines of No. 10 or 14 gauge steel, properly ventilated, 
when lined with 2 inches of red oak or equivalent material, shall be approved. Extra 
precautions that may be taken are: removal of wheels to prevent theft of magazines; 
a second lock on the door as required; a mortice-type lock if obtainable, or when pad- 
locks are used, they should be covered with hoods; or the forming of a barricade around 
the magazine by use of bags filled with sand or dirt. 

2. When the detonator magazines are constructed as provided in Section 2B(c), 
of the Instructions, or as described in your letter, and are securely fastened to a piece 
of drill casing 10 to 15 feet in length, which is inserted in a vertical hole, they may be 
approved as satisfactory. Additional precautions may include addition of a second lock 
as in (1) and the barricading of the box by use of bags filled with sand or dirt. 

3. In regard to shooting trucks: two magazines will be allowed on the same truck, 
one to be used for explosives and the other fo. ‘etonators, when properly constructed 
and lined. They shall be kept locked at all times excepting when in use. A partition 
placed in one magazine to make separate compartments will not be acceptable. Other- 
wise, trucks shall conform to all the other regulations for transportation as provided 
in the Instructions. 

4. Where companies do not at the present time possess portable metal magazines, . 
they may submit to this office, or any district Bureau office, specifications of the 
construction details of the type of magazine that they wish to build, and the office will 
advise whether or not it will meet the necessary specifications. Thought should be 
given to constructing the magazine along the following lines: 

5. A portable magazine constructed of 2 inches of tongue and grooved oak lumber, 
2 inches of packed cotton fibre, and 1 inch of matched flooring on the outside, will be 
acceptable for this type of work. Further recommendations may include the covering 
of the magazine with No. 26 gauge iron, if obtainable. Under certain conditions, if the 
iron cannot be obtained, asbestos or similar fire-resisting material might be accepted 
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as covering for this wood. In all cases it will be necessary that the magazines be in- 
spected by a Bureau of Mines representative before it can be considered as approved. 


I wish to express the appreciation of the Bureau and myself for the cooperation 
you have given us in this work. 
Very truly yours, 
(Signed) G. M. Kintz 
Supervising Engineer 
District G 


APPOINTMENT OF REGIONAL EDITORS 


Section V-B-4 of the Constitution and By-Laws of the Society provides that the 
Editor may appoint associate, regional, and special editors to assist him in performing 
his editorial duties. The Editor wishes to announce the appointment of the following 
Regional Editors: 


Dr. L. L. NETTLETON, Gulf Research and Development Company, Pittsburgh, 
Pennsylvania, for the EASTERN AREA; 

Dr. W. M. Rust, Jr., Humble Oil and Refining Company, Houston, Texas, for 
the GutF Coast AREA; 

Pror. Dart WANTLAND, Colorado School of Mines, Golden, Colorado, for the 
Rocky Mountain AREA; and 

Dr. H. W. WasHBurn, Consolidated Engineering Corporation, Pasadena, Cali- 
fornia, for the Paciric Coast AREA. 


CHANGE OF ADDRESS OF BUSINESS MANAGER 


Mr. J. F. Gate, Business Manager of the Society, has changed his address to 
P. O. Box 1925, Washington, D. C. All communications other than those per- 
taining to editorial matters should continue to be sent to Mr. Gallie, at his new 
address. 


MEMBERSHIP APPLICATIONS APPROVED FOR PUBLICATION 


The Executive Committee has approved for publication the names of the follow- 
ing candidates for membership in the Society. This publication does not constitute an 
election, but places the names before the membership at large. If any member has in- 
formation bearing on the qualifications of these nominees, he should send it to the 
Secretary within thirty days. (Names of sponsors are placed beneath the name of each 
nominee.) 

ACTIVE 

Roland Leroy Clifton 

John A. Gillin, C. G. McBurney, W. B. Hogg 
Dallas Robert Davis 

J. H. Deming, George D. Mitchell, Jr., Roy L. Lay 
Carter Holt Gregory 

W. Hafner, H. R. Thornburgh, Harry L. Thomsen 
Claude Macy Hathaway 

Charles C. Williams, C. G. McBurney, John A. Gillin 
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Oscar Branche Jackson 

F. Goldstone, M. O. Gibson, Claude Wendell Horton 
Harold Hunter Jones 

T. I. Harkins, R. S. Jackson, Joseph L. Adler 
Clarence Henry Leathers 

L. F. Athy, Anton H. Erickson, E. J. Handley 
Daniel Linehan 

M. M. Slotnick, Derry H. Gardner, Joseph Lynch 
Phillip Francis Martyn 

F. F. Reynolds, J. P. Black, J. P. Schumacher 
Mark Howard McKinsey 

F. Goldstone, L. K. Mower, Claude Wendell Horton 
Linden Jay Richards 

Donald Beeth, L. F. Athy, H. R. Prescott 
Theodore Rozsa 

F. Goldstone, L. K. Mower, L. E. Deacon 
David Sheffet 

Joseph A. Sharpe, Daniel Silverman, M. M. Slotnick 
Joseph Carl Waterman 
W. Hafner, H. R. Thornburgh, Harry L. Thomsen 


ASSOCIATE 
Guy O. Buckner, Jr. 

T. I. Harkins, A. L. Smith, R. S. Jackson 
Donald Crary 

D. F. Broussard, T. I. Harkins, R. S. Jackson 
James Greene Jackson, Jr. 

Cecil E. Reel, C. Hewitt Dix, Lawrence K. Morris 
Raymond Morris Low 

T. I. Harkins, A. L. Smith, R. S. Jackson 
Walter Obier Martin 

Walter S. Olson, Louis A. Scholl, Jr., C. C. Zimmerman 
Leroy Robert Meitzen 

T. I. Harkins, A. L. Smith, R. S. Jackson 
Edward Alfs Pielemeier 

Barthold W. Sorge, C. Hewitt Dix, Herbert Hoover, Jr. 
John Harold Ware 

T. I. Harkins, A. L. Smith, R. S. Jackson 
Roscoe C. Wilber, Jr. 
T. I. Harkins, S. P. Weatherby, R. S. Jackson 


TRANSFER TO ACTIVE 
Felix Lemuel Yeates 



































PERSONAL ITEMS 


J. ANGELO SMITH, on leave from the Humble Oii and Refining Company, is serv- 
ing in the Radiation Laboratory of the Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts. 


Car L. Bryan, formerly with the Gulf Research and Development Company, 
has been commissioned a Lieutenant in the U. S. Naval Reserve. Mail will be for- 
warded to him from 2405 Isabella Avenue, Houston, Texas. 


Leo Horvitz has organized the Horvirz RESEARCH LABORATORIES, located at 
3217 Milam Street, Houston, Texas. 


SHERWIN F. KELLEy reports that he has recently been engaged in geophysical 
studies by gravitational, electrical, and magnetic methods on zinc deposits in Ten- 
nessee, graphite in Pennsylvania, and iron in New Jersey. 


Lt. CLEVE C. Nasu, JR., is stationed with Company C, 27th Signal Training 
Battalion, Camp Crowder, Missouri. 


FLoyp WIIL.ts, JR., Research Engineer with Geophysical Research Corporation. 
should be addressed at 2507 N. Boston Place, Tulsa, Oklahoma. 


J. P. Woops, formerly with the Shell Oil Co., Inc., is Research Associate in the 
Radio Research Laboratory at Harvard University, Cambridge, Massachusetts. 


Mark MItsTEt, consulting geophysicist, has opened an office in New York City 
at 489 Fifth Avenue. Although he is now stationed in New York, he retains an office 
at 155 Macpherson Avenue, Toronto, Canada. 


W. B. HeEroy, President of the Pilgrim Exploration Company, has been named 
to succeed E. DeGolyer as director of conservation in the Office of the Petroleum 
Coordinator. He may be reached at the Cosmos Club in Washington, D. C. 


R. B. Moran, JR., is now acting as Research Engineer in the geophysical labora- 
tory of the Superior Oil Company. His address is 3328 Martha Circle, Pasadena 
California. 


N. A. HASKELL is now located at 9 Everett Street, Cambridge, Massachusetts. 


W. D. Mownce, on leave from the Humble Oil and Refining Company, is serving 
in the Naval Ordnance Laboratory in Washington, D. C. His residence is 1921 Kalo- 
rama Road, N.W. . 


Major DREXLER DANA may be addressed c/o Mrs. Dana in Newton, Illinois. 


A. T. SCHWENNESEN has his office at 2702 Niels Esperson Building, Houston, 
Texas. 


G. W. HosTETLER has been advanced to the rank of Captain in the Signal Corps. 
He is now located at 1403 South Main Street, Goshen, Indiana. 


O. C. Otson has completed a six months’ contract with the U.S. Engineers as an 
Associate Engineer and may now be addressed at P. O. Box 349, Brady, Texas. 
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; Captain G. D. LAMBERT, Battery “F,” 63rd CA (AA), is now stationed in Seattle, 
i Washington. 


A. E. McKay, formerly with the Atlantic Refining Company, is a Lieutenant in 
the Fighter Command School at Orlando, Florida. 


J. E. Jett, Gulf Coast supervisor for the Geotechnical Corporation, may be 
reached at 408 Sue Street, Houston, Texas. 


Copies of a new 36-page booklet, entitled ““Du Pont Propucts For SEISMIC 
PROSPECTORS,” may be obtained by addressing E. I. du Pout de Nemours and Com- 
pany, Room 2506-R, Nemours Building, Wilmington, Delaware, while the limited 
supply lasts. The booklet contains instructive material on the best methods of priming 
cartridges, loading charges and firing shots, and lists precautions which should be fol- 
lowed in handling, storing and using explosives. 


Lt. E. L. FETZER should be addressed c/o Mrs. Fetzer, 401 LaSalle Street, Nava- 
sota, Texas. 


Lt. J. P. Faris, Jr., has been transferred from Randolph Field to B. O. Q. Elling- 
ton Field, Texas. 





R. H. TALitMan is serving with the International Petroleum Company, Apartado 
803, Guayaquil, Ecuador. 














Noyes D. Smita, JR., formerly with the Shell Oil Co., Inc., is serving as physicist 
in the Naval Ordnance Laboratory, Washington, D. C. His mailing address is 2100 
34th Street, S.E., Washington, D. C. 







G. C. Howarp, of the Carter Oil Company, has been transferred to Ferriday, 
Louisiana, with the postal address Box 548. 












James C. TEMPLETON, Managing Director of the Geophysical Prospecting Co., 
Ltd., of London, i:ngland, as returned to England from South America and may 
again be addressed at 6 Neville Court, Abbey Road, London, N. W. 8, England. 











M. L. BENKE has left the Texas Company and is serving as a First Lieutenant in 
the U. S. Army Signal Corps. He may be addressed at 802 Waverly Avenue, San 
Antonio, Texas. 








LEo J. PETERS, of the Gulf Research and Development Company, has been trans- 
ferred from Tulsa, Oklahoma, to Pittsburgh, Pennsylvania, with the mailing address 
Box 2038. 






O. F. VAN BEVEREN has left Bahrein Island in the Persian Gulf, where he was 
serving with the California Arabian Standard Oil Company, and is returning to his 
home at 4361 Colfax Avenue, North Hollywood, California. 












Captain P. M. HonneE Lt of the Signal Corps is serving in the Department on 
Chemistry and Electricity of the United States Military Academy, West Point, New 
York. 


GEorGE B. SomERs is on travel for the Bureau of Ordnance of the Navy Depart- 
ment. His forwarding address is 1004 Rivermet, Fort Wayne, Indiana. 
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G. E. Waite should be addressed at 200 St. Paul’s Road N., Hempstead, New 
York. 


L. F. FiscHER has accepted a position as Chief Computer on Party No. 2 of the 
General Geophysical Company. Mail should be sent him at 833 Fairview Street, 
Jackson, Mississippi. 


Marion A. ARTHUR, on leave from the Humble Oil and Refining Company, is 
serving in the Department of Terrestrial Magnetism, Carnegie Institution of Washing- 
ton, Washington, D. C. 


Louris C. PAKISER, JR., recently elected a student member while in attendance 
at the Colorado School of Mines, has accepted employment as Attached Helper with 
the Carter Oil Company, P. O. Box 801, Tulsa, Oklahoma. 


CLYDE H. Witson is acting as a Mine Operator, engaged in producing and shipping 
tungsten ores. His mailing address is Gold Hill, Utah. 


Harry R. IMLE has been advanced to the rank of Captain of Field Artillery and 
is serving on the Staff and Faculty of the Field Artillery School, Communication De- 
partment, Fort Sill, Oklahoma. 


J. W. Tuomas, Jr., Observer with the Stanolind Oil and Gas Company, should 
be addressed at Box 652, Mt. Pleasant, Texas. 


BuLteETIN No. 15 of “The Experiment Station Soil Mechanics Bulletin” issued by 
the U. S. Waterways Experiment Station in Vicksburg, Mississippi, contains a German- 
English vocabulary of soils terms compiled by H. B. Edwards. 


MELvIn C. TERRY is on leave from the Humble Oil and Refining Company and is 
serving in the Naval Ordnance Laboratory in Washington, D. C. 


J. W. Louse, is serving as Engineer in the Applied Physics Laboratory of Johns 
Hopkins University, Silver Springs, Maryland. 


W. B. WILeEy, of the Gulf Research and Development Company, is located at 158 
Ridgewood Drive, Metairie, New Orleans, Louisiana. 


Former CoRPORAL CHARLES C. BaTEs of Battery A, 2nd F. A. Obsn. Bn., is now 
an Aviation Cadet, located at 109124 Wellworth, Village Station, Los Angeles, Cali- 
fornia. 


Bitty E. Ricuarps, of the United Geophysical Company, is located at P. O. 
Box 946, Perryton, Texas. 


Dwicat E. Warp, formerly of the Carter Oil Company, has accepted foreign em- 
ployment with the International Petroleum Company, Box 803, Guayaquil, Ecuador. 


Major O. T. Hatiipay, of the Signal Corps, 8th Corps Area, may be addressed 
at Fort Sam Houston, Texas. 


MICHAEL ALLON and R. G. DE Goop may be addressed c/o Gulf Research and 
Development Company, Pauls Valley, Oklahoma. 
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Major Epwarp Lipson has advised the Society that his ac” ss, until further 
notice, will be 2900 Argyle Drive, Monticello Park, Alexandria, Virginia. 


S. A. MartIn is located at 14 Skahan Road, Belmont, Massachusetts. 


Mavrice SKLAR, of the Shell Oil Co., Inc., has been transferred from Groesbeck, 
Texas, to Wiggins, Mississippi, where his mailing address will be Box 36. 


L. L. ANTES, on leave from the Humble Oil and Refining Company, is serving in 
Engineering Science and Management Defense Training, U. S. Office of Education, 
University of Texas, Austin, Texas. 


CAPTAIN JOHN A. PRANGLIN, Commanding Officer, Co. B., 833 Engr. Bn., may be 
addressed c/o Judge John L. Pranglin, Pearsall, Texas. 


RayMonD A. PETERSON is acting as Group Supervisor of the Columbia University 
Division of National Defense Research. His mailing address is 120 Pierce Avenue, 
East Greenwich, Rhode Island. 


Francis A. ROBERTS, of the Carter Oil Company, has been transferred from 
Hazelhurst, Mississippi, to P. O. Box 622, Norman, Oklahoma. 


E. I. Heap, of the Gulf Research and Development Company, has been trans- 
ferred from Lubbock to Fort Stockton, Texas, with the mailing address Box 881. 


Howarp Itren, of the Stanolind Oil and Gas Company, may be reached at Box 
652, Mt. Pleasant, Texas. 


E. W. JoHNson should be addressed c/o Cities Service Oil Company, Geological 
Department, Pratt, Kansas. 


J. M. KixsBorn, of the Gulf Research and Development Company, should be 
addressed General Delivery, Hastings, Nebraska. 


J. F. GALuiE, Business Manager of the Society of Exploration Geophysicists, has 
joined the procurement division of the Naval Ordnance Laboratory in Washington, 
D. C. His new address is P. O. Box 1925, Washington, D. C. 
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